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ABSTRACT 

The effects of dark matter annihilation on the evolution of intergalactic medium 
(IGM) in the early Universe can be more important if dark matter structure is more 
concentrated. Ultracompact Minihaloes (UCMHs), which formed by dark matter ac- 
cretion onto primordial black holes (PBHs) or initial dark matter overdensity produced 
by the primordial density perturbation, provide a new type of compact dark matter 
structure to ionize and heat the IGM after matter-radiation equality z cq , which is 
much earlier than the formation of the first cosmological dark halo structure and later 
first stars. We show that dark matter annihilation density contributed by UCMHs can 
totally dominated over the homogenous dark matter annihilation background even 
for a tiny UCMH fraction /tjcmh = ^ucMH(Zcq)/£W > 1(T 15 (1 + z) 2 (m x c 2 /100 
GeV)~ 2 / 3 with a standard thermal relic dark matter annihilation cross section, and 
provide a new gamma-ray background in the early Universe. UCMH annihilation be- 
comes important to the IGM evolution approximately for /ucmh > 10~ 6 (m x c 2 /100 
GeV). The IGM ionization fraction x- lon and gas temperature T m can be increased from 
the recombination residual x- lon ~ 1CU 4 and adiabatically cooling T m oc (1 + z) 2 in the 
absence of energy injection, to the highest value of x- lon ~ 0.1 and T m ~ 5000 K at 
z > 10 for the upper bound UCMH abundance constrained by the cosmic microwave 
background optical depth. 

A small fraction of UCMHs arc seeded by PBHs. The X-ray emission from gas 
accretion onto PBHs may totally dominated over dark matter annihilation and become 
the main cosmic ionization source for a PBH abundance /pbh = ^pbh/^dm S> 10~ n 
(10~ 12 ) with the PBH mass M PBH ~ 10 _6 M Q (10 2 M Q ). However, the constraints 
of gas accretion rate and X-ray absorption by the baryon accumulation within the 
UCMHs and accretion feedback show that X-ray emission can only be a promising 
source much later than UCMH annihilation at z < z m <C 1000, where z m depends on 
the masses of PBHs, their host UCMHs, and the dark matter particles. Also, UCMH 
radiation including both annihilation and X-ray emission can significantly suppress the 
low mass first baryonic structure formation. The effects of UCMHs radiation on the 
baryonic structure evolution are quite small for the gas temperature after virialization, 
but more significant to enhance the gas chemical quantities such as the ionization 
fraction and molecular hydrogen abundance in the baryonic objects. 

Key words: intergalactic medium — dark matter — cosmology: theory — early 
Universe — galaxies: structure 



1 INTRODUCTION 

Ultracompact Minihaloes (UCMHs) are primordial dark 
matter structures which formed by dark matter accreting 
onto primordial black holes (PBHs) after matter-radiation 
equality z cq ~ 3100, or direct collapsed onto an initial 
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dark matter overdensity produced by small density per- 
turbation befo re z ea , e.g., in several Universe phase tr an- 
sition epochs (|Mack et al.ll2007l ; iRicotti fc Gouldll2009h . If 
the density perturbation in the early Universe exceed a crit- 
ical value S c — [5p/p) c ~ 1/3, this region becomes grav- 
itationally unstable and directly collap se t o form a PBH 
(|Hawkind Il97ll : ICarr fc Hawkind 1 1974 : see iKhlopovi |2010| 
for a review and references therein). PBHs which form with 
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a sufficient high mass > 10 16 gram do not evaporate but be- 
gin to grow by accreting the surrounding dark matter and 
form a compact dark matter halo, which will grow by two 
ord ers of magnitude in mass during the matter dominated 
era (|Mack et al.ll2007l ) . These haloes are so-called Ultracom- 
pact Minihaloes (UCMHs), or say Primordially-Laid Ultra- 
compact Minihaloes (PLUMs). On the other hand, small 
density perturbation in the early Universe 10 -3 < S < 8 C 
will form a compact dark matter overdensity instead of a 
PBH. Such an overdense cloud can also seed the formation 
of UCMHs llRicotti fc Gouldll2009l ; IScott fc Sivertssonll2009l; 
iJosan fc Greenll2010h . Note that the initial density pertur- 
bations from inflation were just 5 ~ 1CT 4 - 10~ 5 , it is pro- 
posed UCMHs are far more viable to form by accreting onto 
dark matter overdensity, which requires a much lower per- 
turbation threshold than PBHs. Also, the UCMHs seeded by 
primordial overd ensities have a different profile with those 
seeded by PBHs l|Bertschingerlll985l ; [Mack et al.ll2007l ). 

UCMHs have been recently proposed as a new type 
of no n-baryonic massive com pact gravitational object (MA- 
CHO; HicottiXGould 2009) as well as gamma-ray and neu- 
trino source (|Scott fc Sivertssonll2009l '). UCMHs could pro- 
duce a microlensing lightcurve which can be distinguished 
from that of a "point-like" object such as a star or brown 
dwarf, thus become a promising new target for microlensing 
searches. Moreover, the abundance of UCMHs can be con- 
strained by the observation of the Milky Way gamma-ray 
flux and the extragalactic gamma-ray background, although 
this constraint is still very uncertain based on today's data 
(ILacki fc Beacomll2010l ;|j osan fc Gr een 2010; Saito fc Shirail 
l201ll ). Since we kn ow the growth of an isolate UCMH as a 
function of redshift (|Mack et al.ll2007l ) , we can natively trace 
the fraction of UCMHs back to very high redshift without 
considering mergers and tidal destruction. Until now, most 
of the works on UCMHs focus on the properties of the nearby 
UCMHs at z < 1. Another important question which has 
barely been discussed is that, what are the consequences of 
UCMH radiation at very high redshift, since UCMHs are 
the "remnants" originally from the early Universe? As the 
sources of heating and ionization before the first structure, 
stars and galaxies, sufficient UCMHs may play important 
roles to change the chemical and thermal history of the early 
Universe. Our main purpose in this paper is to investigate 
the impacts of UCMH emission on the intergalactic medium 
(IGM) in the Universe reionization era, and the following 
first baryon structure formation and evolution. 

The process of reionization of all hydrogen atoms in 
the I GM would have b e en completed a t redshift z ~ 6 
(e.g., iBecker et al.l 120011 ; iFan et all |2002| ) . However, much 
earlier ionizati on at z > 6 is implied by the WMAP ob - 
servation (e.g.. iDunklev et al.ll2009l ; iKomatsu et al.l f2009). 
It is commonly suggested that the possible contributions 
to the high redshift reionization between approximately 
6 < z < 20 are the first baryonic objects to produce sig- 
nificant ultraviol et light, early (Pop HI and Pop II) stars , 
and old quasars (iBarkana fc Loebl[2007t IWise fc Abelll2008l ; 
IVolonteri fc Gnedinll2009l ; lMeiksinll2009h . However, it is still 
unclear whether quasars and first stars were sufficiently ef- 
ficiency to reionize the universe. Dark matter, on the other 
hand, is suggested to be the exotic source of ionization and 
heating at high redshift due to its self-annihilation or decay. 
It is usually proposed that weakly interacting massive parti- 



cles (WIMPs) provide a compelling solution to identify the 
dark matter component. The mass of the dark matter parti- 
cles m x , and the average annihilation cross section (av) are 
the two crucial parameters to affect the ionizing and heat- 
ing processes. Under the thermal relic assumption that the 
cross section {av) ~ 3 x 10 -26 cm 3 s -1 to match the ob- 
served QomIi 2 w 0.110, most previous studies showed that 
the effects of homogenous dark matter background annihi- 
lation or decay on the high-redshift IGM are expected to be 
important only for light dark matter m Y c 2 < 1 GeV or ster- 



ile neutrinos (e.g., [Hansen fc Haimad 2004; iPierpaoli 2004 
i Mapelli fc Ferraral 120051; iBelotskv et al.ll2005l ; IZhang etaT 



ba 



2006; IMapelli et al.ll2006l ; iRipamonti et alj|2007al lbt IChlub; 
2010). The annihilation flux would be enhanced only af- 
ter the formation of the first dark objects z < 60, as 
dark matter become more clumpy (e .; 



Chuzhovl 120081; 



iNataraian fc Schwarj 120081 . 120091 . |2010| ; iBelikov fc Hooper] 
|2009| , |2010| ). However, in our case, as dark matter is more 
concentrated in UCMHs which are significantly denser than 
the homogenous dark matter background, WIMP dark 
matter annihilation within UCMHs may become powerful 
gamma-ray source dominated over the homogenous back- 
ground annihilation, even though UCMHs are very rare. 

A small frac tion of UCMHs are seeded by PBHs 
l[Mack et al.ll2007l ). In this paper we call these UCMHs as 
PBH host UCMHs. Since PBH abu ndance is still uncer- 
tain for a broad range of PBH mass (|Josan fc Green! 120091 ; 
ICarr et alJlioibl ), we only give a qualitative estimate that 
the abundance of PBH host UCMHs should be much less 
than other UCMHs. For PBH host UMCHs, the X-ray emis- 
sion from the accreting baryonic gas flows onto PBHs may 
totally dominated over dark matter annihilation within the 
host UCMH, since the Eddington luminosity is several or- 
ders of magnitude brighter than that of annihilation from 
the host UCMH, and the photoionization cross section for 
hydrogen or helium is much larger than the Klein-Nishina 
or pair production cross section for energetic gamma-rays. 
It is very difficult for a "naked" PBH to re ach a sufficient 
high accretion rate in the IGM environment. (jBarrow fc Silk! 
1979l;ICarr!ll98ll;lGendin et al.lll995l;lMiller fc Ostrikerll200ll ; 
Ricottill2007l ; lRicotti et al.ll2008l ; iMack fc WeslevH2008r ). but 
the situation will be quite different if PBHs are surrounded 
by UCMHs. The accretion rate and X-ray luminosity of 
baryons can change significantly when the effects of a growth 
UCMH is involved (|Ricottill2007j ; iRicotti et al.ll2008l ). How- 
ever, it is possible that the gas is heated and piles up around 
the PBH if the host UCMH is sufficiently massive. Also, 
the accretion feedback such as outflows or radiation pres- 
sure prevent gas from being totally eaten by the PBH im- 
mediately, if the gas accretion rate significantly exceeds the 
Eddington limit. As a consequence, the gas density and tem- 
perature within the UCMH may be significantly higher than 
the cosmic universal gas density, and the X-ray emission is 
totally absorb in the UCMH, but reradiate basically in the 
infrared band. In this paper we will give criteria for X-ray 
emission escaping from the host UCMH to ionize the IGM. 
Also we will compare the importance of X-ray emission from 
PBH host UCMHs and dark matter annihilation from total 
UCMHs in the early Universe, depending on the abundance 
of both total UCMHs and PBH host UCMHs. 

Another topic related to the UCMH radiation is that, 
the formation and evolution history of the first baryonic 
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structure can be changed by UCMH radiation. Previous 
studies showed that the annihilation or decay of the ex- 
tended distributed dark matter in the first structure both 
change the gas temperature and the chemical properties 
such as t he abundance of molecular coolants such as H2 
and HP dBiermann fe Kusenkoll2006l ; IStasielak et al. I l2007l ; 
iRipamonti et alT2007bl ). Higher coolants abundance helps 
to decrease the gas temperature and favors an early collapse 
of the baryon gas inside the halo, but dark matter energy in- 
jection delays this collapsing process. It is still under debate 
whether dark matter annihilation or decay inside the dark 
halo will promote or suppress the first structure formation. 
Nevertheless, it is concluded that either the promotion or 
suppression effect is quite small for most dark matter mod- 
els, as the change of gas temperature in a virialized halo 
for various dark matter models is small. If the first large 
scale dark haloes contain UCMHs, these UCMHs can in- 
ject more annihilation energy into the halo than the first 
dark haloes, and potentially play more important role to 
change the properties of the first haloes than the extended 
distributed dark matter in haloes. Moreover, X-ray emis- 
sion which comes from PBHs also suppress the formation of 
the first baryonic objects. Therefore it is also worthwhile to 
study the effects of UCMH radiation on the first structure 
formation and evolution. 

This paper is organized as follows. In Section [2] we 
calculate the dark matter annihilation luminosity from 
UCMHs, and the X-ray emission from PBH gas accretion. 
We emphasize on the importance of UCMH annihilation 
compared to the homogenous dark matter background an- 
nihilation, and focus on the physical reasons that whether 
and when the X-ray emission from PBHs becomes more im- 
portant than the UCMH annihilation in the early Universe. 
In Section [3] we discuss the gas heating and ionization pro- 
cess by the two types of UCMH radiation from 2 ~ 1000 to 
10, and investigate the impact of UCMH radiation on the 
IGM evolution. Next in Section [4] we show the influences of 
UCMH radiation on the first baryonic structure formation 
and evolution. The main results of this paper are given in 
Sections [3] and [4] In Section [5] we discuss the importance of 
UCMHs in the reionization era, the effects of a single massive 
UCMH in the first baryonic structure, as well as other sec- 
ondary effects. Reader could skip this section and directly go 
to Section[6] which presents the conclusions. In this paper we 
do not consider dark matter decay, which should have similar 
consequences as the annihilation process. Also we fix the an- 
nihilation cross section to be thermal relic {av} — 3 x 10 -26 
cm 3 s _1 , although much larger cross section (av) — 3x 10 -24 
cm 3 s _1 to 10~ 20 cm 3 s _1 is proposed in the hope of ex- 
plaining the reported Galactic cosmic r ay anomalies as the 
results of dark matter annihilation (e.g. , Chang et alJl200"S 
lAbdo et al"1l2009l ; IXharonian et alj|200g| ). A larger cross sec- 
tion with the same dark matter particle mass m x can have 
higher luminosity and more significant influence on ionizing 
and heating the early Universe. Table [T] gives the notation 
and definition of some quantities in this paper. 



2 RADIATION FROM UCMHS 

In this section we discuss two types of energy emission from 
UCMHs in the early Universe: dark matter annihilation, and 



X-ray emission from the accreting baryonic gas onto PBHs. 
As mentioned in Section 1, the second type of emission is 
related to a small fraction of UCMHs, which host PBHs. 
Generally we still call the second type of energy emission as 
UCMH radiation, that is because a PBH is always located in 
the center of its host UCMH and belong to a PBH-UCMH 
system. 

2.1 Dark Matter Annihilation 

Dark matter annihilation lumin osity of nearby UCMHs 
(z = 0) is calculated recently (IScott fc Sivertssonl 120091; 
Lacki fe Beacomll2010l :lj osan fc Greenll2010l ; ISaito fc Shirail 

10 



201 lj) . 



We assume that UCMHs stop growing at 



when the structure formation progressed deeply to pre- 
vent dark matter from further accreting. Now we calcu- 
late the annihilation luminosity as the function of red- 
shift in the early Universe before z ~ 10, and compare 
the result to the homogenous background annihilation. The 
mass of the UC MHs accreted by dark matter rad i al in- 
fall is given by dMack et al l 120071; iRicotti fc Gould] 120091 : 



IScott fc Sivertssonll2009l ; I Josan fc Greenll2010h 



mh(z) — 8m 



1 + z cq 
l + z 



(1) 



where z cq « 3100 is the redshift of matter-radiation equality, 
and 8m is the mass of initial dark matter overdensity. The 
density profile in an UCMH oc r~ a can be written as 



Px{r,z) 



(3 - a)m h (z) 



(2) 



where the factor (3 — a) /47r in equation ((2j) is obtained by 
normalizing the total mass inside the maximum halo extent 
radius Rh as 8m, and Rh is calculated by 



R h (z) « 0.019 pc 



1000 
l + z 



( m h (z) 



1/3 



(3) 



Dark matter annihilation reduces the density in the in- 
ner region of an UC MH, and makes th e density in this region 
to be flat. Following rUllio et all l|2002h the UCMH power-law 
density distribution is truncated at the maximum density 



p(r C ut 



where t 



Pmax 



M(t-t,)' 



-3/2 



(fir, 



(4) 

^^Hq 1 is the age of the Uni- 



verse at a certain redshift z, and t\ ~ 77 kyr is the initial age 
at z oq . Thus the total dark matter annihilation luminosity 
within the UCMH can be calculated as 

2 2 / \ / \ 2 



L 



2nr n x (r)(av)m x c dr 





2ttc 2 



2a 



K-{av) — (l + zY 



where K 



3 V 2q -3- 
xSm(t — ti)~ 3 m x a , 
= (3 - a) (4.66 x 10 8 ) Q " 3 (1 + z c 



(5) 



W3, 



(47r) cgs 

units, n x (r) = p x (r)/m x is the dark matter particle number 
density. The UCMH density pr ofile can change from a steep 
slope a = 3 (Mack et al. 2007)) for the outer part region to 
a = 1.5 ()Bertschingerl 19851 ') for the inner part region, if there 
is a PBH in the center of the UCMH. In particular, radial 
infall onto a central extended overdensity shows a profile 
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Notation and definition of some quantities in this paper 



notation 


definition 


§/Eq. 




977,1. (z\ 


UCMH mass at redshift z 


EH 


eq 




Px ( r ; z ) 


UCMH (density) profile at z 


EH 


eq 


l2l 


Rhiz) 


extent radius of UCMH 


EH 


eq 


3 
1 * u 


m x 


DA^l particle mass 


EH 


eq 


L *J 




DM average anni cross section 


eh 


eq 


m 


J, 


3 n fi i lnm r\T a tMTitrlo TT(;A/TrT 


EH 


eq 




/uCMH 


/iJCMH(^eq) ^TJCMH (-^eq )/^DM 


EH 


eq 


[71 


^ann 


UCM^H anni lum per volume 


EH 


eq 


[To] 


L bkga 


homogenous DM! anni background 


EH 


eq 


[Til 

L__J 


Z a cc 


X-ray rad density from PBH host UCMHs 


EH 


eq 


pL3l 


JPBH 


\ AID T3 HJ /,\ / nAf 


EH 


eq 


[131 


1 13 


Bondi accretion radius 


EH 


eq 


|T9l 


A 


amplification factor of the IGM T m 


EH 


eq 


|T9l 


77 1 


Dimensional accretion rate 


EH 


eq 


[231 


.1 b 


baryonic fraction in a UCMH 


EH 


eq 


[24] 




charac. redshift for gas accretion 


EH 


eq 


[28l 


x- (z) 

■"'ion \^ ) 


reionized baryon fraction at z 


EH 


eq 


[35] 


t(z) 


energy deposition rate per volume at z 


EH 


eq 


[BUI 




photon energy from anni DM 


EH 


eq 


|4T1 


Ex 


charac. energy for X-ray emission 


EH 






T m 


IGM gas temperature 


EH 


eq 


[43] 


fa 2 


molecular hydrogen fraction in IGM gas 


clo i 1 
MJj 


eq 


H7J 


-^halo 


total anni lum from a dark halo 


EH 






^UCMH 


anni lum from UCMHs inside a halo 


EH 








anni lum from extended DM in a halo 


EH 






e loc,iso 


energy deposited by isothermal DM halo anni 


EH 


eq 


Us] 


e loc,UCMH 


energy deposited by UCMHs anni inside a halo 


EH 


eq 


Ha 


^loc,acc 


energy deposited by X-rays inside a halo 


EH 


eq 


EH 


„halo / bkgd 

z m I ^rn 


z m for PBH host UCMHs inside/outside a halo 


EH 








energy deposition rate density by <5-func SED 


EH 


eq 


EH 


factor 


ratio of rad of differently distributed UCMHs 


EH 


eq 


EH 



Table 1. In this table "DM", "anni", "lum", "charac", "rad", "func" are short for dark matter, annihilation, luminosity, characteristic, 
radiation and function. The subscript "acc" is for X-ray emission because X-rays are emitted by gas accretion onto PBHs. 



p oc r , which is more widely use d as the typical density 
profile for most region i n UCMHs (iRicotti fc Gould! 12003: 



- - regit. - - .. 

IScott fc Sivertssonlboogl : IJosan fc Greenll2010T i. Taking a 
9/4, we have 



Lann = 36.3L (ffv) X J 3 171^/^ ( 1 + z) ( 



(6) 



where (av} s — (av}/3 x 10 -26 cm 3 s _1 and m x ,ioo = 
m x c 2 /100 GeV. According to equation the annihilation 
luminosity decreases with the evolution of the Universe, ba- 
sically because the annihilation flats the inner density pro- 
file as showed in equation Q. The left panel of Fig. [1] gives 
UCMH annihilation luminosity with different halo profile 
a and dark matter particle mass m x . Note that the anni- 
hilation luminosity can be much brighter for lighter dark 
matter particles, and a shallower density profile reduces 
iann significantly!]] In the limit case a ~ 1.5 or 3, we have 



1 In Fig.[T]the annihilation luminosities following the density pro- 
file equation J2} are somewhat overestimated for a steep profile 
a ~ 3, because in this case the halo mass within the truncated ra- 
dius r cu t can no more be neglected. Thus the normalization factor 
of the density profile is oc [ln(i?^/r cu t)] _1 , which is different from 
the factor (3 — a) /(An) in equation J2}. However, as we show that 
L alm for a steeper UCMH profile leads to several orders of mag- 
nitude higher L ann than that with a = 2.25, the conclusion that 



i ann oc (av)/m x or L ann to be independent with m x and 
(av) . The abundance of UCMH as the function of redshift 
is still uncertain today. It can be presented by a parameter 



/uCMh(z) = fiuCMH(z)/!2DM, 



(7) 



where J-tjcmh and Qum are the comoving abundances 
of UCMH and total dark matter with SIucmh(z) = 
f2ucMH(zcq)(l + Zcq)/(1 + z) and Q DM (z) = ^dm(z = 0). 
We have /ucmh(z = 0) ~ /ucmh(z«i)(1 + «eq)/(l + 10) ~ 
3 x 10 2 /ucMH(z cq ), which shows that the UCMH mass grows 
by up to two order of magnitude from z cq to z ~ 10. From 
now on we take /ucmh as /uCMH(z eq ) for short to show 
the initial abundance of UCMH at matter-radiation equal- 
ity, and in current stage /ucmh is taken as a parameter for 
simplicity. 

The mean free path of gamma-ray photons from an 
UCMH with energy _B 7 is written as 

1000 \ 3 



> 



3xl0 3 pc 



nA{z)a(E 1 ) UA(z)ar ' \1 + z J 

where ua(z) = ua(1 + z) A is the atomic number density at 
redshift z. The average distance of inter-UCMH is estimated 



a steeper profile gives a brighter annihilation will not be changed 
too much even in the limit case a = 3. More details of the UCMH 
profile are discussed in Section 15.21 



(8) 
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2.2 Gas Accretion onto PBHs 




100 

redshift 



10 1000 100 
redshift 



Figure 1. Left: UCMH dark matter annihilation luminosity with 
the halo profile a = 1.5 (black lines), 2.25 (blue lines), 2.9 
(red lines), and dark matter particle mass m x c 2 =100 GeV (solid 
lines), 1 GeV (dashed lines), 100 MeV (dotted lines). We adopt 
5m = IMq in this figure. Right: ratio of UCMH luminosity to 
homogenous dark matter background annihilation, where we take 
the fraction of UCMH in the total dark matter as /tjcmh = 10~ 4 
at z = 2 oq , and the lines as the same as in left panel. 



rfuCMH = 



1 


1/3 


(6m) 


1/3 


nucMH(z)_ 




_/ucmhPdm(z)_ 





7r>c f _1/3 

' JUCMH, 



(Sm) 



1/3 



1000 

1 + 2 



(9) 



We have the mean free path exceed the inter-UCMH dis- 
tance Aucmh 3> rfuCMH except for an extremely small 
/ucmh < 10~ 12 . Therefore the cosmic UCMH annihilation 
also gives an uniform gamma-ray background radiation field 
as well as that produced by the homogenous dark matter. 
UCMH annihilation luminosity per volume is given by 

Zann = ^7"" /uCMH (-Zcq)pDM ( Z = 0)(1 + z) 3 

dm 

i a „ m-28 -3 -1 / vl/3 -1/3 r /i i \4 

= 1.4x10 ergem s (<rv) s ' m Xil ' 00 /ucMH(l + z) . 

(10) 



On the other hand, the energy injection rate by the self- 
annihilation of the homogenous dark matter background per 
volume is 



Zbkgd = (av) 



2m v 



= 3.2x 10" 43 ergcm- 3 s- 1 ( ( j U ) s m- 1 100 (l + 2) 6 . (11) 

We compare the radiation between UCMH and normal dark 
matter annihilation as 

: 4.5xl0 12 ( ( ™)J 2/3 m 2 / 3 00 (^Tl) /uCMH » (12) 



'bkgd 



If /ucmh > 2.2xlO- 15 ( ( 7i;) 2/3 m x 2 1 / 3 (l+ Z ) 2 , the gamma-ray 
background due to dark matter annihilation is dominated by 
UCMH annihilation. More details depending on the density 
profile a and dark matter m x can be seen in the right panel 
of Fig. [1] 



The abundance of PBHs Opbh as a fraction of total dark 
matter Qdm at z < z cq can be parameterized as /pbh = 
fipBH/^DM. We ignore the PBH growth and take /pbh as 
a constant in the matter-dominated Universe for two rea- 
sons. The first reason is that, as the accretion processes had 
been significantly suppressed before z ~ 10 due to the rela- 
tive motion between PBHs and baryon gas, the PBH growth 
timescale is i gr0 wth ~ ^Saip — 5 x 10 s yr just reaches or is 
longer than the universe age at t(z ~ 10) ~ 5 x 10 8 yr. The 
second reason is that low mass PBHs has lower accretion 
rate while high mass PBHs are inclined to produce outflows, 
which further increases the accretion timescale and makes 
PBH growth to be negligible compared to its host UCMH 
growth. A very simi lar statement to kee p a constant /pbh 
was also proposed in lRicotti et all l|2008l ). Keep in mind the 
PBH abundance /pbh is different from the UCMH initial 
abundance /ucmh at z cq as mentioned in Section 12.11 be- 
cause a large amount of UCMH seeds at z eq should be the 
initial primordial dark matter overdensity but not PBHs. 
According to the density primordial perturbation theory, 
generally we have the relation /pbh "C /ucmh, which will 
be discussed in details in Section \2. 2. II 

Much work has been done to show the effects of radi- 
ation from PBH or early black hole accreti on on the early 
Univ e rse thermal and ionization history (e.g. .[Barrow fc Silk 
19791: iGendin et al.l 1 19951: iMiller fc Ostrikerl IgOOll ; IRicotti 
20071 : IRicotti et al.ll2008l ; lRipamonti et al.ll200Sl ). Our goal in 
this section is to focus on the importance of PBH gas accre- 
tion radiation compared to the overall UCMH dark matter 
annihilation. The X-ray emission form accreting PBHs may 
lead to a very different heating and ionization history of the 
early Universe compared to the dark matter annihilation. 
The X-ray luminosity from an individual PBH with mass 
Mpbh can be written as r^Lsdd = 'i-Kr/GmpMpBu/o'TC ~ 
3.3 x 10 3 77_iLq(M P bh/Mq) with L E dd and = 77/O.I 
being the Eddington luminosity and average radiation ef- 
ficient of all PBHs respectively. This X-ray luminosity is 
much higher than the dark matter annihilation luminosity 
equation ©. Thus the X-ray radiation density in the early 
Universe z > 10 can be written as 



A/pBH 

1.3 x 10 



i?/pbhPdm(z = 0)(1 + z) 3 



erg 



cm 3 s 1 77_i/pbh(1 + z) 



(13) 



Combing equations (|10[) and (|13[) . the ratio between PBH 
accretion luminosity and UCMH dark matter annihilation 
luminosity is 



9<ffw>: 



-1/3 1/3 

m x,ioo 



T] -if PBH 

/ucmh 



10 

l + z 



(14) 



Since the IGM heating rate due to the energy injection by 
PBH X-ray emission or UCMH dark matter annihilation is 
proportional to both the energy injection rate, and the IGM 
cross section for all the interactions suffered by the UCMH 
emitted photons in X-ray band Ex or annihilation emit- 
ted gamma-ray photons E ly the importance of the IGM gas 
heating by X-ray emission and UCMH dark matter annihila- 
tion can be estimated by the ratio Z a ccftot(Sx)/iannO'tot(S 7 ) 
with rjtot labeling the total cross sections in different pho- 
ton energy range. If we roughly take the X-ray and IGM 
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interaction cross section atot(Ex) the Thomson cross sec- 
tion, and the high energy photon interaction atot(Ej) the 
Klein-Nishina cross section (see Section 13.11 for more accu- 
rate calculations), the energy deposition in the IGM due 
to gas accretion e acc and annihilation e ann radiation is esti- 
mated as 



Can 11 



-V3 4/3_ 



V-ifp 



X.100 I f 

V /UCMH 



10 

1 + z 



-(15) 



which gives the first conclusion that the X-ray heating may 
become totally dominated over dark matter annihilation in 
the early Universe if the PBH abundance exceeds a critical 
value as 



> 3.1 x 10' 

/UCMH 



/ \l/3 -4/3 

(av) s ' m x ,ioo 



v 10 



(16) 



Theoretically the value of nfpBH //ucmh includes many 
uncertainties. In general, there are at least three reasons to 
have a low value ?7 /pbh //ucmh <C 1: density perturbation 
scenarios prefer low initially value of /pbh//ucmh; ineffi- 
cient radiation 77 1 is favored by low mass PBHs while 
accretion feedback decreases rj for high mass PBHs or PBHs 
with high mass UCMHs; and X-ray emission from PBHs 
can be trapped inside the surrounding host UCMHs which 
accumulate baryons. 



2.2.1 PBH Abundance 

Either PBHs or UCMH overdensity seeds are produced ba- 
sically by the density perturbation in the very early Uni- 
verse during some special epoches such as inflation or phase 
transitions. The cosmological abundance of UCMHs can be 
estimated by integrating from the overdensity seed thresh- 



old 



10" 



to the PBH formation threshold 5 C 



1/3 

(jRicotti fc Gouldll2009l;IScott fc Sivertssonll2009l) . Similarly, 
the PBH abundance is e stimated by integrating the pertur- 
batio n above 5 C ~ 1/3 (I Green fc Liddlelll997l ; iGreen et all 
Il997l ). Assuming a Gaussian perturbation at a formation 
redshift Zf 3> z aq to produce both PBHs and UCMH seeds, 
the ratio /pbh//ucmh at matter-radiatio n equality can be 
directly traced back to formation time zj l|Carr et al.ll2010l ; 
lKhlopov|[2oioh . As a result, the relative abundance of PBHs 
to UCMH overdensity seeds formed at redshift Zf is written 



/pbh 
/ucmh 



2a(z f Y 



d5 



exp 



0X P [-2^ 

ro- 6 -<5r 



(15 



2a 2 



exp 



18a- 



(17) 



The perturbation variance at zf is r oughly given by cr(zf) ~ 
9.5 x 10- 5 [M h or(«/)/10 66 g] (1 "™ )/4 l|Green fc Liddld 1 1997ft . 
with Mhor(z/) and n being the hori zon mass and mass spec- 
trum index at Zf. Taking n < 1.3 ()Lidsev et al.lfl995h . the 
ratio /pbh //ucmh from a Gaussian perturbation is the func- 
tion of horizon mass as 



fp 



/ucmh 



< exp 



M hol (z f ) 
5.5 x 10 10 e 



(n-l)/2' 



(18) 



which means the value of /pbh //ucmh becomes <C 1 for 
Mhor{zf) 2> 5.5 x 10 10 g, not to mention the fact that 



the masses of dark matter overdensity seeds or PBHs are 
even lower than the horizon mass 5m <C Mhor(zf) and 
ilfpBH <S M hol (zf). Combing equations (fT5")l and (JTSJ, X-ray 
emission from gas accretion hardly becomes the dominated 
heating source in the early Universe, except for low mass 
PBHs AfpBH < M hor (z f ) < 3.7 x 10 18 g in the Gaussian 
perturbation scenario. However, PBHs in this mass range 
should either have disappeared within a Hubble time due to 
the Hawking evaporation, or too small to accrete the IGM 
gas. 

As a result, the initially Gaussian density perturbation 
at a certain epoch is not able to generate sufficient abun- 
dant PBH to dominated over the total UCMH dark matter 
annihilation emission, basically because the large amplitude 
part of a Gaussian distribution is highly suppressed. On the 
other hand, non-Gaussian perturbation may give an even 
lower probability of PBH formation, as the large fluctuation 
can be suppressed in the non-Gaussian d istribution and fur- 
ther decrease the ratio of /pbh//ucmh (|Bullockfc Primackl 
Il997h . 

However, other mechanisms such as different formation 
epoches for UCMHs and PBHs, different early inflationary 
potential, double inflation models, various phase transitions, 
and cosmic string collapse may enhance the high ampli- 
tude perturba tion and increase the PBH abundance (see 
iKhlopovl |2010| and references therein). Also, it is still ar- 
guable if all the 5 > 10~ 3 perturbation could produce dark 
matter overdensity in the radi ation dominant era. For ex- 
ample, iRicotti fc Gouldl (|2009h requires the host UCMHs 
around PBHs t o have similar initia l pertu rbation amplitude 
as PBHs, while [Scott fc SivertssorJ (|2009h has less strict re- 
quirement as 5 > 10 -3 to form the initial dark matter over- 
density. There are more physical uncertainties to estimate 
the abundance of PBHs and UCMHs produced by other 
mechanisms than a simple Gaussian distribution assump- 
tion. Therefore we still take /pbh as a free parameter which 
satisfies /pbh <C /ucmh to describe the relative abundances 
between PBHs and UCMHs. 



2.2.2 Inefficient Radiation 

Another effect to constrain the X-ray luminosity density in 
the early Universe by the gas accretion onto PBHs is the 
low radiation efficiency due to low accretion rate onto low 
mass PBHs, or the significant radiative feedback, thermal 
outflow and suppressed accretion rate due to accretion onto 
high mass PBHs or PBHs with high mass host UCMHs. 

In principle the mass distribution of PBHs is broad 
enough to cover the range from the Pla nck mass ~ 10~ 5 
g to thousands of solar mass 10 s M Q (e.g.. lCarr et"al]|2010h . 
As mentioned in the above Section 12.2.11 if PBHs are 
formed from the Gaussian perturbation with the variation 
a oc M~ {n ~ 1)/4, and the index n > 1, low mass PBHs 
should be more abundant because of the higher density per- 
turbation variance a for lower mass M. Also, phase tran- 
sition models give PBH mass or UCMH seed less than 1 
Mq l|Scott fc Sivertssoni r2009'). On the other hand, keep in 
mind in the IGM environment a "naked" PBH without a 
host UCMH can never reach the Eddington accretion rate 
Af£>p d H = L Edd /c 2 ~ 1.4 x lO 17 (AfpB H /M ) g S -\ unless its 
mass is M PBH > 360 M Q [1000/(1 + z)f /2 . The surrounding 
host UCMH increases the accretion rate if the PBH mass 
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is A'/pbh > 100M Q l|Ricotti et al.ll2008l , their Fig. 4). Note 
that an ideal case is r\ ~min{0.1m, 1} after the accretion be- 
come super-Eddington, while the typical accretion efficiency 
for quasars or microquasars disk is 77 ~ 0.15. For low mass 
PBHs with m < 1 the radiation effici ency is est imated as 
77 ~ 0.01m 2 for the spherical case (|Shapirolll973al lbh. which 
gives much lower efficiency than the high accretion rate that 
V < 0.1. 

If high mass PBHs (WOMq < A/pbh < 1O 5 M ) suc- 
cessfully form with an appreciable abundance compared to 
the lo w mass PBHs, a s discussed by some previous au- 
thors (|Mack et al 1 l2007l ; ISaito" et al. 1 120081 ; iFrampton eFal] 
120101 ), or the host UCMH seeds are more massive than 
the PBHs Sm > M PBH (|Ricotti fc Gould! 120091 , more de- 
tails see Section f2. 2. 3p . the Bondi accretion rates onto these 
PBHs with their host UCHMs can significa nt exceed the 
Eddi ngton limit after some critic redshifts (|Ricotti et al.l 
2008). However, a spherical super-Eddington accretion 
is generally unstable and inclined to drive high mass 
loss rate with thermal outflows (e.g., ISmith fc Owockil 
2006). Recent simulations show that radiative feedback 
may become important to reduce or even quench t he 
accretion process pe riodically l|Milosavlievic et all l2009al lbl; 
I Park fc Ricotti feoill ). Also, the thermal heating by the out- 
flow energy or radiative feedback will increase the tem- 
perature of the gas around PBHs and decrease the Bondi 
radius and accretion rate onto PBHs. Besides the spher- 
ical accretion case, the falling gas angular momentum 
will become important for m 2> 1, and form an ac- 
cretion disk around PBHs. However, the physics of the 
super-Eddington accretion disks is still not clearly known. 
Various types of super-Eddington accretion disk models 
have been proposed, such as the o ptically-thick advec- 
tion dominated accr etion flow fADAF. iNaravan fc Yi|[l994l . 
Naravan et al.l fl998l) , the ad iabatic inflow-outflow (ADIO, 
Blan dford fc Be gclman 1999), the convection -dominated ac- 



cretion flow (CD AF. INaravan et al 



doughnut s" torus Jib ramowicz et al 



.2000]) . the "polish 
Il978l ) and the thick 



slim disk (|Abramowicz et al] |T988). In most cases the super- 
Eddington accretion disk advects most of its heating en- 
ergy inward into the black hole without emission, and 
has a lo w radiation efficiency r? for high accre t ion r ate 
r\ < 1 (lAbramowicz et all 1 19881 ; INaravan et al.l Il998l or 
[Ab ramowicz fc FragilelboTll for a review). 

In a brief summary, low average radiation efficiency n in 
equations (|13|) to (|15[) are favored because of the low accre- 
tion rate onto the low mass PBHs, and radiative or viscous 
feedback and outflows of accretion onto high mass PBHs or 
PBHs with high mass host UCMHs, which also leads to a 
low value of >)/pbh//ucmh and suppress the importance of 
PBH X-ray radiation from PBHs compared to the overall 
UCMH dark matter annihilation. From now on we consider 
the X-ray emission is mainly contributed by the accreting 
PBHs with m > 1. 

2.2.3 Radiation Trapping in Host UCMHs 

Some previous works discussed that the accretion flow 
around PBHs is Compton thin in most cases, since in the 
sub-Eddington accretion case the spherical flow is transpar- 
ent near the PBH, while in the super-Eddington accretion 
case the accretion flow are inclined to form an accretion 



disk (e.g.. lRicotti et al]|2008l ). However, sufficient high mass 
UCMHs can accrete and thermalize baryons from the am- 
bient IGM, even there are no PBHs in the center of these 
UCMHs. As the gravity potential at the outer edge of the 
host UCMH is mainly contributed by the UCMH mass, but 
the accretion onto the center PBHs is according to the PBH 
mass, the accretion rate into the host UCMHs is not neces- 
sarily equal to the accretion rate onto the center PBHs. In 
other words, baryons can be firstly accumulated and viri- 
alized inside the host UCMH during the accretion from 
the IGM to the inner UCMH region, followed by a sec- 
ondary accretion onto the center PBH and feedback (out- 
flow) from the accreting PBHs. Based on this consideration, 
the baryons inside the UMCH can be divided into two com- 
ponents: the piled up baryons inside the UCMH, and the 
accretion spherical flow or disk around the center PBH. Al- 
though the optical depth of the accretion gas or disk, which 
is mainly contributed by the depth around the inner horizon 
region r ~ i?sch is transparent to X-ray photons, the X-ray 
emission can still be trapped and absorbed by the piled up 
baryons inside the host UCMH, and reradiate photons with 
much longer wavelength into the outer IGM environment. 
Quantitative analysis is given as follows. Part of the treat- 
ment is similar to an analogy discussion on the dark matter 
structure formation and ba ryons filling process |Hoefl et al.l 
120061 ; lOkamoto et"ai1l2008l ) . 

If UCMH dark matter annihilation does not change the 
IGM temperature evolution, the IGM temperature is ap- 
proximately coupled with the cosmic microwave background 
(CMB) temperature before the decoupling time Zdcc ~ 100, 
and the IGM sound speed before z& cc is c s ~ 5.7 km 
s _1 (^j§) 1//2 ■ In general, the UCMH annihilation heating 
and PBH emission without trapping increases the IGM 
temperature. We introduce an amplification factor A that 
T m = -4Tcmb at z > Zdoc , where T m and Tcmb are the tem- 
perature of the IGM and CMB respectively, and A depends 
on the UMCH profile and annihilation properties, as we will 
calculate in Section [3] The sound speed c s oc T 1 ^ 2 becomes 
c 3 ~ 5.7 km s~ 1 A 1 / 2 (y^) 1 ^ 2 , an d the Bondi accretion ra- 
dius (i.e., the accretion sonic sphere) of a PBH-UCMH sys- 
tem at z > Zdcc is 



Gm h 400 pc 1 
TB ~ ~cT ~ (ITT) 2 



5m 

Mo 



(19) 



Equation (|19[) is derived under the assumption that the 
Bondi radius is larger than the UCMH size rs > Rh- Fur- 
thermore, if tb > 2Rh, the virial temperature of the host 
UCMH T vir ~ (wf) [ Gm R h {Z) ] is greater than the tem- 
perature of the ambient IGM gas T v i r > T m . According to 
the general virial theorem, the thermal pressure of the gas 
due to virialized heating is weak compared to gravity of the 
UCMH. In this case we consider the IGM baryons should fall 
into t he UCMH unimpeded, regardless of th e center PBH 
mass (Hocft ct al. 2006; lOkamoto et~aT1l2008l ). 
The criterion rs > 2Rh at 2 > z& cc gives 



1^ ) > 1600^ 2 ( i±£ 
M ) V 1000 



(20) 



Note that higher IGM temperature around UCMH, i.e., 
higher A gives a higher minimum UCMH mass to attract 
baryons. Similar result can be derived for the case after de- 
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coupling z < Zdec, where the IGM gas temperature decou- 
pled with the CMB temperature and dropped adiabatically 
as T a( j oc (1 + z) 2 without any heating sources. We still take 
the factor A > 1 to measure the IGM temperature increase 
due to annihilation T m = AT a d- Then using the criterion 
tb > 2Rh, we find that baryons fall into UCMHs unim- 
peded at z < Zdcc if 

5/2 

(21) 



(Mucmh - MpbhX* - ti) ~ m h (z)f b 



(24) 



> 16M 



3/2 ( l+_g 

100 



As a result, if the UCMH initial overdensity seed is 8m > 
16OO.4 3/2 M for z dcc < z < 1000, or 5m > 160^ 3/2 M Q for 
z < Zdcc, the IGM gas can always fill into the UCMH no 
matter it includes a PBH or not. Otherwise for a lower 8m, 
the critical redshift z c below which the UCMH accrete is 
(1 + z c ) < 0.63.4~ 3/2 (<5m/M©) for z > z dcc and (1 + z c ) < 
13A- 3/5 (Sm/M Q ) 2/5 for z < z dec . If the UCMH hosts a 
PBH in the center, baryons are still able to piles up and 
thermalized in the host UCMH due to the gas virialization. 

The lower bound of gas accretion rate into the UCMH 
can be estimated as 



Mucmh = 47rr 2 3 c s/ 9g as (z) > \-kR\v / ;(Rh)m b n h (z) 

~ 8.2xl0 16 gcm- 3 (l + ,) 1 / 2 (|g), (22) 

where Vff(Rh) is the free fall velocity at Rh- If all the gas 
into the UCMH is totally accreted onto the center PBH, the 
dimensionless accretion rate of the PBH is m = M g /Mj^\& 
is 

rh> l 8 .0['i+£V / 7^)>l, (23) 



1000 



Mp 



with M™ being the Eddington limit accretion rate onto 
the central PBH. Note that the ideal rh can be even higher 
if the initial host UCMH is more massive than the cente r 
PBH 8m » AfpBH as discussed in iRicotti fc Gould! (2009). 
However, the real accretion rate should be less than the 
value in equation (|23[) for two reasons. First, the baryons 
can be heated and virialized during the accretion process 
in the UCMH, and has a temperature ~ T v i r warmer than 
the IGM T m to increase the gas pressure and decrease the 
accretion rate onto the PBH. And the gas temperature is 
further increased 3> Tvir near the PBH due to PBH emis- 
sion and ionization. Also, super-Eddington accretion disks 
are also inclined to drive outflows. The positive Bernoulli 
parameter over most of the ADAFs due to t he small radia- 
tion loss may trigger strong outflows or jets (|Naravan fc Yil 
Il994h . and produce an ADIO s in which outflow carries awa y 
most of flow mass and energy IjBlandford fc B cgclman 1999). 
Also, CDAFs may produce a "convective envelope" with no 
accretion onto the black hole (|Naravan et al-lfeoOOT ). In gen- 
eral accretion disks with super-Eddington accretion rate are 
inevitably accompanied by outflows and winds, which signif- 
icantly decrease the final accretion rate onto the black hole. 
In the PBH case, these outflows should be injected back to 
the host UCMH environment. 

The upper bound of the baryonic fraction in the UCMH 
is the universal fraction Q,\,/Q. m . However, as the UCMH 
grows following equation ([T}, we adopt a more conservative 
method to estimate the lower bound of baryon fraction /j 
inside the UCMH. We estimate the baryonic fraction in the 
UCMH fb (the mass ratio between gas and dark matter) as 



where we take Mucmh S> Mpbh, i.e., most of the accreted 
gas into the UCMH is piled up without being immediately 
eaten by the PBH. Combining equations (122 p and (|24|l . we 
have the lower bound of baryon fraction to be > 7.6 x 
10 -3 , which is a constant independent of the redshift. 

The optical depth of the piled up gas in the UCMH due 
to Compton scattering is 



/im p 



(25) 



Since ft from equation (|24|l is a constant, and UCMH growth 
does not change the steep region pdm oc r~ a with a = 
9/4 but only increase Rh and flats the region r < r cu t (see 
equation [4]), we take the baryon fraction to be uniformly 
distributed in the UCMH, both in the steep and flat region. 
Therefore the column density of the baryon gas inside the 
UCMH depends on the UCMH profile, which depends on 
the dark matter properties {(av), m x ) given by equation Q. 
Actually the baryon profile can be steeper in the flat region 
of the halo r < r cut since the dark matter annihilation flats 
the inner halo profile, thus gives an even larger optical depth. 
Furthermore, we consider the baryon gas is ionized x e ~ 1 
inside the UCMH, at least in the flat region r < r cu t- We 
check that the Stromgren radius of the PBH emission rs 
satisfies r cut < rs < Rh, as the heated gas near the PBH can 
reach a temperature as high as the Compton temperature 
~ 10 keV in the ionized region. The hot ionized gas around 
the PBH produces a small sonic sphere in the dense baryon 
region near the PBH, decreases the accretion rate onto the 
PBH, giving Mucmh S> Mpbh as mentioned in equation 
((24}. Note that there should be two distinct sonic spheres, 
the sphere for the host UCMH outside Rh, and that for 
the center PBH i nside the UCMH. This scenario is similar 
with lWang et all |2006l ) that an accreting BH has two Bondi 
spheres, a smaller inner sphere in the hot gas region and 
a larger one in the outer cooler region. Therefore we take 
x e ~ 1 in equation (|25[) . The optical depth is written as 

I °~Tfb 



> 10 



p x (r)dr 



( 6m\ fl + zY' 5/9 1-13/9 / f\ 

ImoJ (TuOoJ <iooW S 7 - (26) 

Hereafter we take {av) s — 1. Combining equations (|20[) and 
(|26|l . we conclude that before the decoupling z > Zdcc the 
gas is always Compton thick to the X-ray emission from the 
center PBH accretion when the host UCMH itself accretes 
baryons. After the decoupling z < z dc c the redshift range 
that X-ray emission escapes is 



(1 + «)< 13 A- 3 '* < 7m~J« A^ 

(1 + *)< 62 (^) " 2/5 m-%1 (If) > 7m-% 6 A*/ 4 



From equation (|27p there is a maximum z m in the case of 
r B > 2R h that 



(r B >2R h ) 



1/3 j-3/10 



(28) 



if z > Zm, X-ray photons will be totally trapped. Note that 
z m insensitively decreases with the increasing of the IGM 
temperature factor A. A larger optical depth due to a deeper 
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baryon profile at r < r cut gives an even lower z m . Also, the 
range of Sm applied in equation (|27[) is 



0MM e A 3/2 <Sm< 80M Q m~% 3 . 



(29) 



In other words, in the r B > 2Rh case, no X-rays can escape 



the host UCMH if Sm > 9,QM G) m 



-5/3 



© m x,ioo- 



More massive PBHs 



are easier to reach Eddington accretion rate, but more dif- 
ficult to produce a transparent baryon environment in the 
host UCMHs. 

On the other hand, if r B < 2R h (i.e., T vir (R h ) < T m ), 
most part of the UCMH gravity potential well is not deep 
enough to compress the gas and overcome the pressure 
barrier of the gas virialization heating. In this case the 
UCMH itself cannot accrete and thermalize baryons except 
for the region in the radius r' B where satisfies [Gmh(r < 
r' B ) J (2r' B )] — c 2 . At z > Zd cc , the critical radius r' B for a 
pure UCMH profile p x oc r~ a is 

l/(c-2) 



r' B = Rh 



Gm h 
2dR h 



(30) 



where we apply a = 9/4 from equation ([5]). The part of 
UCMH inside r' B can accrete and heat baryons. Similar to 
equations (|22p to (|26p . the accretion rate into the region 
r < r' B in the unit of Eddington accretion rate of the center 
PBH at z > z^c is 



2ivr B c s mbnb(z) 



M PBH 

iu Edd 



3.9 x 10~ 1X A 



-15/2 



100 

1 + z 



9/2 



8m 

M© 



dm 



4) 



where the factor 2-n is due to the suppressed accretion at 
r > r' B in the UCMH, thus the baryon density is half of the 
ambient gas density. Assuming 8m — Mpbh, equation (|31[) 
shows that only high mass PBHs (A/pbh > IOOMq) are able 
to produce super-Eddington accretion if there is no accre- 
tio n feedback. This is basically consistent with the results 
in iRicotti et al.l (120081 ) . However, we should mention two 
things. The first thing is that, if 8m > A/pbh, the ideal accre- 
tion rate in equation (|31[) also increases. The second thing, 
which is similar to the analysis below equation (|23p is that, 
the real accretion rate onto the PBH is lower than the ideal 
m due to higher gas temperature and accretion feedback. As 
a result, we find that baryons can be accumulated and virial- 
ized inside the UCMH region r < r' B with a baryon fraction 
approximately as f b ~ 10~ 3 ^~ 9/2 (l+z)~ 3 (<5m/M Q ) 3 . Using 
the condition m > 1 in equation (|3ip at z > z^cc for a suf- 
ficient radiation efficiency, the baryon optical depth inside 
radius r' R is 



i.i x io- 2 A- V2 {i + zy LA/a m a <i 



13/6 5/9 



*' 100 ^ \ Sm ) 



8m \ 10/3 

5/6 



LOnW^P ( l -±l) >l, (32) 



100 



where we also take the baryon profile inside the UCMH 
is proportional to the dark matter profile for simplicity, 
and Xe ~ 1. According to equation (|32p . we consider the 
UCMH is optically thick to the X-ray emission in the case 
of r B < 2Rh and z > Zdcc, unless 8m 3> Mpbh or dark 
matter particle mass m x ,ioo <IC 1. 

After decoupling z < Zdcc, we find that baryons can be 
accumulated and virialized inside a radius 



r' B » 0.012 pc A 



30 y / 8m 
with a upper bound z m that 

1/2 



z£ B<2Rh) ^32A 



3/8 / Sm \ _ 5/12 _ 



(33) 



(34) 



below which (z < z m ) the gas is Compton thin inside ra- 
dius r' B . Higher ratio <5m/MpBH S> 1 or lighter dark matter 
wi Xi ioo «C 1 increases z m , which also decreases slightly if 
annihilation effect is included to heat the IGM gas (A > 1). 

Not only the baryons inside UCMHs can trap X-ray 
photons, but also the outflows driven by PBH accretion feed- 
back also absorb the X-ray emission as well. As mentioned 
before, the super-Eddington accretion rate onto the PBH 
is unstable to trigger strong outflows in both spherical and 
disk cases. An optically thick "outflow envelope", both in 
the polar and equatorial region around the PBHs, forms to 
cover the PBH and totally or mo stly absorb X-ray emission 
form the inner accretion flows (llgumenshchev et al.1 [20031; 
Kohri et all 120051 ; iPoutanen et al.1 120071 ; lAbolmasov et al.1 
2009). In this case, the Compton heating in the outflow re- 
gion should be important to increase the gas pressure and 
temperature, balance the gravity well, reemit thermalized 
photons from outflows, and reg ular and the accretion rate 
onto the PBH (|Wang et al.ll2006r ). It is likely to have a steady 
state or periodically changing outflow envelope covering the 
whole PBH, but the details are still an open question which 
is beyond the purpose of this paper. What we want to show is 
that, even though the accretion disk itself around PBH is op- 
tical thin to X-ray radiation, the X-ray emission can be still 
absorbed in the outflow envelope due to the accretion feed- 
back and disk instability in the super-Eddington case, not 
to mention the optical-thick or geometry thick disks which 
absorb X-ray emission by themselves. 

We give a summary of Section 12.21 Equation (|13p com- 
putes the X-ray radiation density due to baryon gas accre- 
tion onto PBHs. The ratio ?7/pbh//ucmh is parameterized 
in this paper. The much lower probability of PBH formation 
compared to the UCMH formation, and inefficient radiation 
due to low mass PBHs (Sm <C IOOMq ) or accretion feedback 
from high mass PBHs or PBHs with massive UCMHs (Sm S> 
Mpbh) give ??/pbh//ucmh -C 1. Moreover, we simply in- 
troduce a critical redshift z m below which (z < z m ) X-ray 
emission from the super-Eddington accretion PBHs (m ^> 1) 
become important to heat and ionize the early Universe. 
Hotter IGM heated by other energy source s (e.g. , annihi- 
lation) slightly decreases z m . IRicotti et al.1 (|2008T > showed 
that the accretion becomes m > 1 at z m ~ 20 (100) for 
Mpbh = 10 2 (300) M , and m is always m 3> 1 (m <C 1) 
for Mpbh > 1O 3 M (M PB h < 10M©). However, as we dis- 
cussed in Section \2. 2. 31 the stage of super-Eddington accre- 
tion heating the IGM can be delayed, because X-ray photons 
are trapped inside the total or inner region of the UCMHs 
due to the accumulation and virialization of the accretion 
gas. Outflows can also (partly) absorb X-rays. We take the 
critical redshift z m as showed in equations Q28p an d (|34p . 
for 2 < z m X-ray photons from the super-Eddington ac- 
cretion PBHs could escape their host UCMHs. If the PBH 
abundance is much less than that of UCMH, but still satis- 
fies equation (|16p . X-ray emission should dominate over the 
early dark matter annihilation at z < z m . 
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3 REIONIZATION AND HEATING OF THE 
IGM 

3.1 Basic Equations 

The evolution of baryon ionizati on fraction Xj on (z) is given 
by the differential equation (e.g.. ICirelli et al1l2009l ) 



dL 



n A (l + z 



3 dXjonjz) 

dt 



I{z)-R{z), 



(35) 



where I(z) and R(z) are the ionization and recombination 
rates per volume respectively, n A is t he atomic number den- 
sity to day. We use the rate R{z) from lNataraian fc Schwar3 
(2008). The ionization rate per volume due to dark matter 
annihilation or X-ray emission from gas accretion is given 
by 



/(*) 



E x dn(z) 

dE^ t ^-P{E 1 ,z)N io ^{E 1 ), 



(36) 



E x (l + z)/(l + z eq ), the differen- 



J x — "°X L 

ted photon, E cq 
tial term dn(z)/dE- l is the p hoton spectral num ber den- 
sity at redshift z . We follow | Cirelli et al.l |2009| (see also 
iBelikov fc Hooperl 120091 ; iNataraian fc Schward 120081 , 120091 ) 
to calculate the probability of primary ionizations per sec- 
ond P(E 1 , z), and the number of final ionizations that gen- 
erated by a single photon of energy i? 7 produces Ni on (Ej). 
Note that JVion(-E 7 ) is propo rtional to the ionization fac- 
tor r7i on (iri on ) » (1 — Xio n )/3 (|Shull fc van Steenberd fl985l ; 
IChen fc Kamionkowski 1120041 ). which means approximately 
1/3 emitted energy goes into the reionization of atoms if 

Xion < I- 

First of all, we consider the ionization is due to UCMH 
dark matter annihilation. The spectral number density is 
obtained as 



rfw ann (z) 
dE~, 



dz 



/ Cdt dU nn (z') ( 1 + Z 

dz 1 E x dEy(z') \l + z' 



exp(-r), (37) 



where Z ann is the annihilation luminosity as mentioned in 
Section [2TT1 Ey(z') = E 7 (l + z')/(l + z). The optical depth 



r = I dz"^n A (l + z") 3 a tot (E^ 



(38) 



with Eyi = Ey(l + z")/(l + z). The total cross section 
cr(tot) for the DM annihilation photon to interact with 
electrons in the IGM mainly includ es the Klein-Nishina 
cross section for Compton scattering l|Rvbicki fc Lightmanl 
120041 ') and the photonionization cross section for H and He 
as <7H+H e jZdziarski fc Svenssonlll98l^ ). Pair production on 
matter becomes important for m x > 1 GeV. CMB photons 
also contribute to the total cross section for m x > 10 TeV, 
which can be neglected in our cases. 

The total energy deposition per second per volume at 
redshift 2 is given by 



«(*) = 



dE- 



dn a 



dE 



(z) 

^n A (l + 2) 3 o- to t(£ 7 )£; 7 . 



If we take the monochromatic dark matter annihilation 
emission for simplicity, i.e., the photons produced by dark 
matter annihilation are the rest energy of the dark matter 
particle m x c 2 , the photon flux spectral density then can be 
calculated by the <5-function 



dE> 



-CO ^ knn(z')5 (El f - E x ) 



Thus we have the energy deposition e(z) as 

1 + 2 



e(z) 



E * dE„, 4 cdt , 

n A (l + Z) -p-Lnn(Z ) 

-Ecq * 



l + z' 



where z' satisfies 
z' = ^L(l + z)-l, 



(40) 



r a tot (£ 7 ),(41) 



(42) 



and r is calculated from z' to 2. Keep in mind in the 
above formula (I41|l the dark matter annihilation products 
are simplified as the gamma-ray photons with sole energy 
E x = m x c 2 . More realistic annihilation spectrum is model- 
depended. F or example, dn/dE-y can be chosen as fo l lowin g 
the model in lBergstrom et all (| 19981 ) and lFeng et~al] l|200ll) . 

For the X-ray photons from an accreting PBH, equa- 
tion (|41|l will still be available if we choose a to t as the 
X-ray total cross section a to t ~ °"h+h + or, and E x 
with the X-ray ch aracteristic energy Ex as Ex — 3 keV 
(M P bh/M )" 1/4 l|Salvaterra et all 120051 ). As the real ac- 
creting PBH sp ectral energy distribution i s very model- 
depended (e.g., Shakura fc Sunvaevl [ l 973; S az onov et al.l 
|2004 ISalvaterra et al.ll2005l ; iRipamonti et al.|[20ol ). in the 
very first calculation we simplify the X-ray emission as the 
single- frequency emission at a characteristic energy Ex, 
which mostly can be considered as the peaked energy in 
the real spectral energy distribution. We will also discuss 
the more realistic PBH spectral energy distribution in the 
discussion section [B~3l 

Now we list the heating and cooling processes in the 
IGM. The heating of IGM by UCMH annihilation or X-ray 
emission can be written as 



dT m 
~~dt 



Vhent (Xio 



3fcs n A (l + z) 3 (l + / Hc + x ia 



e(z), (43) 



where the heating fraction ?7heat(:Eion) which shows the por- 
tion of energy e(z) into heating IGM is adopted as 77h cat = 
C(l - (1 - x^)) b with C = 0.9971 a = 0.2663 and 
b = 1.3163 dShull fc van Steenberd [l985h . We can approx- 
imated take the He fraction in the IGM as /h g — 0.073. 
Moreover, CMB photons can be treated as another heating 
source for the IGM if the IGM gas is colder than the CMB 
(T m < Tcmb), otherwise the IGM gas would transfer energy 
into the CMB environment. The coupling between IGM gas 
and the CMB photons can be importan t when the differ- 
ence between T m and Tcmb is significan t (|Wevnmannlll965l ; 
iTegmark etHII 19971 ; ISeager et al1 [2000) 



~~dt 



3 2cMB£ion(2cMB — T n 



(44) 



where the coupling rate coefficient fe c 



5.0 x 10' 



Other IGM cooling terms are dominated by the adia- 
batic cooling during the expansion of the universe as 



(39) ( 

dz 



2T m 
I + 2' 



(45) 



for low temperature. Note that the IGM temperature will 
decrease independently as T m oc (I + 2) 2 for a pure adiabatic 
cooling process. Furthermore, for sufficient high temperature 
~ 10 4 K, the molecular hydrogen H2 cooling will also be 
important. This cooling term can be calculated as 
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-TT- ] = Ah, (1 - Xion - 2/h 2 )/h, [riA(l + z) 

v dt /h, 



(46) 



where we adopt the s pecific coolin g coefficient A h , from 
iHollenbach fc McKeel l| 19791 ) and IVoshida et ail (|2006l ). 
We neglect other chemical cooling processes such as 
Bremsstrahlung, helium line cooling, H2 line cooling and 
hydrogen three-body reaction . HP cooling is imp ortant for 
T < 200 K and low density (|Yoshida et al.ll2006h . but the 
gas adiabatic cooling will be dominated in this case. As we 
mainly pay attention to the evolution of the ionized fraction 
Xi on and IGM temperature T m , we only include the evolu- 
tion of hydrogen (H, H~, H + , H2) and electron gas (e~) 
as the main species for ionization. More detailed simula- 
tion including other species and cooling processes is beyond 
our purpose of this paper. The evolution of the H2 fraction 
is ado pted from the semi-analytic model in iTegmark et al.l 
<| 19971 ) 



= kmTlA(l + z) 3 (l - 



2/ H2 )a;io 



(47) 



dfn 2 
dt 

where we follow ITegmark et al.1 (|l997l) and iGalli fc Palial 
(1998) to calculate the reaction coefficient fc m . 



3.2 Solutions of IGM Evolution 

Since no direct evidence related to the UCMH radiation has 
been confirmed until now, the UCMH abundance is still un- 
certain. In this paper we take the UCMH fraction /ucmh as 
a free parameter. In Fig. [2] we show the ionization fraction 
Xi on (z) and the IGM temperature T m for /ucmh = 10~ 4 
and 10 -6 , which correspond to today's expected abundance 
~ 1% and 10 _ respectively. Moreover general discussion on 
the UCMH abundance will be given later. The initial Xi on at 
z — 1 000 is adopted as 0.01, and T* m as the CMB tempera- 
ture llGalli fc Pallalfl99Sl ; iRipamonti ll2007l ; iRipamonti et all 
l2007al lbl). The first basic conclusion which is similar to the 
previous works is that, lighter dark matter particles or 
higher UCMH abundance give larger Xi on and higher T m . 
An extreme bright UCMH annihilation background (e.g., 



m x c < 1 GeV for /ucmh 



10~ 



/ucmh — 10 ), even gives a monotonically increased Xi on 
and T m > 10 4 K without a standard reionization epoch in 
the early Universe. We compare the UCMHs annihilation 
results with the homogenous background annihilation, note 
that the homogenous dark matter annihilation background 
only produce noticeable effects for light dark matter parti- 
cles m x c 2 < 1 GeV or sterile neutrinos. UCMHs, which pro- 
vide a new dominated dark matter annihilation gamma-ray 
background as showed in Section al! play a more important 
role to ionize and heat the early Universe. 

Furthermore, as the cosmological Jeans mass m.j can 
be taken as an indicator of the IGM structure evolution, 
the left panel of Fig. [3] gives the evolution of Jeans mass in 
the cosmic UCMH annihilation background. The Jeans mass 
mj oc Tm 2 p _1//2 should be a constant if the gas tempera- 
ture is always equal to the CMB temperature T m = Tcmb- 
In this paper we call this constant as "CMB mass". In 
the left panel of Fig. [3] mj with various T m is generally 
normalized in the unit of "CMB mass". The remaining 
Thomson scattering optical depth contributed by UCMH 
annihilation is showed in the right panel of Fig. [3] For 



6 < z < 30, the remaining CMB optical depth is estimated 
as St ~ 0.046±0.016 by the WMAP five-year measurement^ 
We assume a linear increase of Xion from z — 10 to the full 
ionization time z — 6, thus the upper bound contribution of 
UCMH annihilation to the measurable CMB optical depth is 
St ~ 0.028 ± 0.016 < 0.044. Based on this consideration, in 
our examples only one extreme case that m x c 2 = 100 MeV 
with /ucmh = 10~ 4 is ruled out by the CMB remaining op- 
tical depth in Fig. UCMH annihilation can significantly 
increase the Thomson optical depth in the early Universe 
2 3> 100 up to St ~ 0.5 without stringent constraints by the 
CMB optical depth measurement at z < 30. After the last- 
scattering epoch the ionization fraction Xi on can change from 
Xion ~ 10~ 4 (without dark matter annihilation) to a upper 
bound Xi on ~ 0.1 (e.g, m x c 2 = 1 GeV and /ucmh = 10 -4 ). 
Also, the IGM can be heated from a temperature T m of adi- 
abatically cooling T m oc (1 + 2) 2 in the absence of heating 
source to the upper bound T m > 10 3 K with a sufficient 
amount of heating contributed by UCMH annihilation. An 
much higher Jeans mass than the "CMB mass" by ~ 2 — 3 
orders of magnitude increase can be obtained due to the hot- 
ter IGM temperature. Therefore, we can natively estimate 
that the formation of small baryonic objects can be strongly 
suppressed, although more investigations need to be carried 
out in Section [4] 

In general, we find the impact of UCMH annihilation 
on the IGM evolution can be empirically estimated by the 



factor m. 



x,ioo/ucmh, 



while the threshold of UCMH abun- 
dance to affect the IGM evolution is approximately given 
by m~j 00 /ucMH > 10 -6 , with an upper bound constrained 
by the CMB optical depth at late times z < 30 as Xion ~ 0.1 
and T m ~ 5000 K at m~^ 100 fvcMH ~ 10~ 2 . The CMB opti- 
cal depth enhancement at early times z > 30 can be more 
dramatic than the late times due to the higher annihila- 
tion luminosity in early redshift (equation [5]), which is dif- 
ferent from the PBH radiation which has higher lu minos- 
ity at late times (Section [2331 iRicotti et afll2008l ). Fur- 
ther phenomenological constraints should be made by CMB 
polarization anisotropies, which is left for a future investi- 
gation. Keep in mind another cannel to concentrate dark 
matter rather than primordial density perturbation is the 
formation of the first dark objects, which should affect the 
IGM evolution much later (z < 100) than UCMHs. There- 
fore UCMH annihilation has a definitely much earlier and 
more important impact on the IGM evolution from the last 
scattering to the structure formation time. 

So far we give the results only for UCMH dark mat- 
ter annihilation. Whether the X-ray emission from the PBH 
host UCMHs will significantly change the above results 
mainly depends on the fraction of PBH /pbh, the aver- 
age inflow radiation efficiency r\ and the critical redshift 
z m as given in Section 12.21 In our paper we combine the 
factor ?7/pbh as one. Remember the results in Section 12.21 
when X-rays from the center PBH region successfully passes 



2 The WMAP five-year measurements give th e CMB Thom- 
son s cattering optical depth r ~ 0.084 ± 0.016 (Komatsu ct al. 
|200S|), which is mostly due to ionization at lat e times z < 30 
llRicotti et~aLll2008l ; iNataraian fc Schwarj|2010l ). If we subtract 
the optical depth contributed by the totally ionized gas t(z < 
6) = 0.038, the re maining depth is St ~ 0.046 ± 0.016 < 0.062 
JCirelli et alJl200Sh . 
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Figure 2. Effects of UCMH dark matter annihilation on the 
IGM evolution. Upper two panels correspond to the case of 
/uCMH = 10 -4 and lower to /tjcmh = 10 — 6 . The thick (black) 
lines from the top down show the results for UCMH annihilation 
with m x c 2 = 100 MeV, 1 GeV, 10 GeV, 100 GeV and 500 GeV, 
while the thin (blue) lines give the results of homogenous dark 
matter background annihilation background for m x c 2 =100 MeV, 
1 GeV and 10 GeV. 




1000 100 10 10 100 1000 



redshift 



/pbh below the value in equation (|16fl • We will give a lower 
limit of /pbh, above which X-rays have obvious impact on 
the IGM evolution at z < z m . In the following calculation 
we assume equation (|16|l is always satisfied, and do not dis- 
tinguish the redshift which divide the UCMH radiation into 
annihilation dominated or X-ray radiation dominated from 
the redshift which gives a transparent baryon environment 
in UCMHs, but simply use one parameter z m . 

In Fig.|4]we take z m and ?i/pbh as parameters with the 
characteristic emission frequency as Ex = 1 keV, 10 keV 
and 100 keV, which correspond to the typical PBH mass 
as 10 2 A/ Q , lO~ 2 Af and 10~ 6 A/ Q respectively. Higher z m 
means higher ratio <5m/A/pBH or lighter dark matter parti- 
cles. Only X-ray emission as the energy source is calculated 
in this figure 0. As showed in Fig. [4] the final properties of 
the IGM at z ~ 10 with same 77/pbh and Ex are more or less 
closed to each other regardless the value of z m , which means 
lower z m gives more dramatic thermal and chemical change 
at z < z m - According to equation (|34[1 . lower z m corresponds 
to lighter m x , which also increase the UCMH annihilation. 
On the other hand, Xi on and T m vary for more than three or- 
ders of magnitude from E x ~ 100 keV (10" 6 M Q ) to E x ~ 1 
keV (1O 2 M0) with the same 77 /pbh, which means massive 
PBH favors the IGM ionization. On the other hand, the X- 
ray radiation effect can be neglected when ti/pbh < 10 -11 
for Ex ~ 100 keV, but a smaller limit r\ /pbh < 10 -12 is ap- 
plied for Ex ~ 1 keV. Below the lower limit the PBHs are 
not expected to have any promising effects on reionization. 
The estimate in Fig. [5] shows that no strict constraints are 
made for ?7/ucmh < 10~ r by the remaining CMB depth St, 
which allows a dramatically increased Jeans mass due to the 
hot IGM gas ~ 10 4 K. 

As a result, X-ray emission from PBHs gives a more 
promising impact on the IGM evolution if 77/pBH S> 10~ n 
(10~ 12 ) for A'/pbh ~ 10~ 6 M Q (10 2 M Q ), or empirically to 
say, 77 /pbh > 1.8 x 10 _12 (A//A/ Q )~ 1/8 . As we assume 
/ucmh < 10 -4 , we expect the UCMH dark matter anni- 
hilation only played its role on the IGM evolution at very 
high redshift z m < z < 1000, but X-ray emission changes 
the Universe reionization history dramatically at relatively 
lower redshift z < z m . Considering 77 ~ 0.1, the upper bound 
value /pbh < 10~ 6 is two orders of m agnitude higher than 
the upper PBH abundance < 10" s in IRicotti et all (|2008T) 
for Mpbh > 10 3 Mq, but muc h lower than the lo w mass 
PBH abundance constraint in IRicotti et all i|200Sl ). How- 



Figure 3. Left: ratio of Jeans mass to the "CMB mass", which is 
the Jeans mass for the gas temperature always being equal to the 
CMB temperature. Right: evolution of the remaining Thomson 
scattering optical depth St with gray belt showing the WMAP 
5-year lc of remaining CMB optical depth at 10 < z < 30. The 
thick (black) lines are for /tjcmh = 10 — 4 and thin (blue) for 
/ucmh = 10~ 6 . The lines with same color from the top down are 
for m x c 2 = 100 MeV to 500 GeV as in Fig. [2] 

through the transparent baryon medium in the host UCMH 
at z < z m , the much brighter X-ray luminosity and much 
larger interacting cross section a to t{Ex) compared to the 
annihilation luminosity and atot(E^) usually guarantees the 
X-ray emission to be dominated over the UCMH annihila- 
tion (equations I15[) , except for a much lower PBH fraction 



3 The more realistic case is that both z m and Ex are the func- 
tions of the mass of PBH and its host UCMH, thus in principle 
neither 2 m nor Ex are single free parameters. The real X-ray radi- 
ation background due to the PBH accretion happens at a certain 
z m as showed in equations J28H and 1134 1 , and change its spec- 
tral e nergy distribution as a function of redshift (Ripamonti et al. 
2008). Technically it is able to calculate the X-ray emission vari- 
ation by introducing the initial mass function of PBHs and the 
host UCMHs. However, both of these two initial mass functions 
are poorly known. More elaborate models only make our model 
more complicated and uncertain. What we focus on in our cal- 
culations is the key differences between X-ray radiation and the 
much earlier occurred UCMH annihilation, so we only use the 
characteristic energy Ex, and the ratio »7/pbh//uCMH to show 
the importance of the X-ray radiation. More discussion of PBH 
spectral energy distribution can be seen in Section 15.31 
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Figure 5. Ratio of Jeans mass to the "CMB mass" (left panels) 
and evolution of Thomson scattering depth (right panel) for the 
starting injection redshift z m = 100 (upper panels) and 20 (low 
panels). The lines are as the same in Fig. [4] 
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Figure 4. Effects of gas accretion around PBHs on the IGM 
ionization and temperature evolution. The upper panels are for 
the X-ray energy injection starting at z m = 100, middle panels for 
Zm, = 50 and lower for z m = 20. The PBH fraction rj/pjn = 10~ 7 
(black lines), 10~ 8 (red lines), 10 -9 (green lines) and 10 — 10 (blue 
lines). Characteristic emission energy are 100 keV (solid lines), 10 
keV (dashed lines) and 1 keV (dotted lines). 

ever, t he massive PBH abundance constraint in lRicotti et al.l 
(2008) is made by the Compton y-parameter estimate at 
Zrcc < z < z cq without local UCMH trapping, while we con- 
sider the two-step accretion first by the host UCMHs and 
then by the center PBHs as mentioned in Section [2.2.31 X- 
rays can just locally heat the accreted gas inside UCMHs 
but not the entire cosmic gas at high redshift. 

The abundance of IGM molecular hydrogen /h 2 in var- 
ious UCMH radiation models is showed in Fig. [6] in this 
section. We see that when the UCMH energy injection can 
be neglected, this fraction goes back to / h, ~ 10~ 6 , which is 
consistent with the standard result (e.g.. lGalli fe Pallalll998l 
and references therein). The upper bound of enhanced /h 2 
is /h 2 ~ 10 -3 , either due to the allowed UCMH dark matter 
annihilation constrained by the CMB optical depth, or the 
X-ray emission at z m < 100. 



4 FIRST STRUCTURES 

In the hierarchical cold dark matter (CDM) scenario, the 
first cosmological objects are dark matter haloes, which 
are formed by gravitatio nal instability from the scale-free 
density fluctuations (e.g.. lGreen et aill2005l ; [Piemand et al.l 




redshift redshift 

Figure 6. H2 abundance evolution for UCMH annihilation (left 
panel) and gas accretion onto PBHs (right panel). The lines in 
the left are the same as in Fig [3] The lines in the right panel are 
the same as in Fig. [4] but thick and thin are for z m = 100 and 
z m = 20 respectively. 

120051 ; lYoshidall2009h . The formation of the first baryonic 
objects depends on the detailed gas dynamical processes. 
The first baryonic objects can successfully collapse and form 
inside dark haloes when its cooling timescale for dissipat- 
ing t he kinetic ener g y is m uch shorter than the Hubble 
time. iTegmark et al.l (ll997T ) showed that the formation of 
the baryonic structure cruciall y depends on the abundance 
of the molecular hydrogen /h 2 ■ iBiermann fc Kusenkd (|2006l ) 
considered that the photons emitted by dark matter annihi- 
lation or decay inside the halo can boost the production of 
H2, and may favor the formation of the first structure. On 
the other hand, a different conclusion that the dark matter 
annihilation or decay can slightl y delay the first baryoni c 
structure formation was given bv lRipamonti et al.l (|2007bl ). 
As they discussed, the higher central density in the baryonic 
cloud without dark matter energy injection could compen- 
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sate the lower abundance of H2 and still lead to the fastest 
cooling. Nevertheless, no matter promotion or suppression 
caused by halo extended dark matter annihilation or decay, 
such effects are pretty small. 

UCMHs can also be captured by the first dark mat- 
ter objects, in this case the UCMH radiation can be much 
brighter than that from the extended large dark matter halo, 
even the fraction of UCMHs in the dark halo is tiny. In 
this section we fucus on the dark halo structures with mass 
~ 1O 6 M (or 10 7 M @ in the PB H heating case), as they favor 
the l ater first star formation rtBroom et all 120091 ; lYoshidal 
2009). For a typical 10 6 Mq halo in the early Universe, 
~ 10 20 (10 18 ) or ~ 10 s (10 6 ) UCMHs within this halo can be 
expected for an initial UCMH fraction /ucmh ~ 10 -4 (10 -6 ) 
if the seed of UCMHs are generated in the electroweak 
or Q CD phase transitions respectively (|Scott fc Sivertssonl 
2009). In these cases the UCMH emission can also provide 
a new type of radiation background in the dark halo. How- 
ever, the uniform UCMH distribution treatment will break 
down if the number of massive UCMHs in a halo is less 
than ~ 10. This happens for the massive UCMH case that 
/ucmhMdm /mh < 10. In this section we study the effects of 
UCMH radiation on the first structure formation and evo- 
lution with /ucmh A/dm 2> fnhi the case of only-several- 
luminous-UCMH will be discussed separately in Section [5]4] 
as a supplementary. As the virial temperature of ~ 1Q 6 Mq 
haloes is less than the threshold for atomic hydrogen line 
cooling, these haloes are often referred as "minihalo" in 
literatures. However, for clearly, in this paper we call the 
first dark matter structure as (cosmological) dark matter 
haloes or dark haloes, which should not be confused with 
the UCMHs. 

4.1 Dark Matter Annihilation 

The profiles of the cosmological dark matter haloes are cho- 
sen before our calculation. The equations for the halo profile 
are listed in the Appendix. Before the formation of the first 
stars, the energy injection inside a large dark matter halo 
mainly contributed by the local emission from the UCMHs 
in the halo, the local annihilation or decay of the extended 
dark matter within the halo, and the outside radiation back- 
ground which injects into the halo. We check the total energy 
produced by dark matter annihilation with a dark halo as 
/-halo = /'UCmh + Lext, with /ucmh and £ cxt being the an- 
nihilation luminosity from UCMHs and the extended dark 
matter in this halo. The typical ratio LuCMH/£ext is demon- 
strated in Fig. [7] where we adopt an isothermal dark halo 
model and /ucmh = 10 -6 . We find both /ucmh and L cxt are 
proportional to the total halo mass A/dm , so /ucmh / L eK t is 
independent to A/dm- In Fig. 0Lucmh is mostly dominated 
in the halo, except for large z v ir with light dark matter parti- 
cles^. g., < 10 GeV for /ucmh = 10~ 6 and z v i r = 100, much 
lighter dark matter particles is required for more abundant 
UCMHs or smaller z v - a ). Therefore, similar to the IGM envi- 
ronment, the energy injection mechanism inside a dark halo 
which contains UCMHs can be very different from the no- 
UCMH case. We focus on the ionization and heating inside 
the dark matter halo. Also, we mention that the results for 
iuCMH/iext in NFW haloes are very similar to the isother- 
mal haloes in Fig. [7] 

The simple criterion of baryonic matter filling in a dark 
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Figure T. Ratio between Lucmh and Lext 

for rn x c 2 =1 GeV 

(solid linos), 10 GeV (dashed lines), 100 GeV (dotted lines) and 
isothermal extended dark matter halo. Thick and thin lines cor- 
respond to the virial redshift z v - lr = 20 and 100 respectively. 
The UCMH fraction is adopted as /ucmh = 10 — 6 . The case 
of //ucmh = Zcxt is marked as the horizon thin dashed line. 
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Figure 8. The minimal dark halo mass for T m = T v i r with the 
IGM heated by UCMHs. The shaded area is for the case that 
the minimal halo mass for gas collapsing becomes Mdm < 6.1 X 
1O 3 M0 due to a lower ambient gas temperature compared to the 
CMB temperature. The lines from m x c 2 = 100 MeV to 500 GeV 
as in Fig. [3] 

halo is approximately taken as the IGM gas temperature be- 
ing cooler than the virial temperature of the halo T m < T v i r , 
otherwise the gas pressure prevents the gas from collaps- 
ing. We do not include the tem perature cooling slope crite- 
rion as m iTegmark et"afl (|l997h to further study the baryon 
cooling and collapsing in the dark halo. Fig. [8] gives the 
dark matter halo mass for the critical case T m = T v i v with 
UCMH annihilation heating the ambient IGM gas. Gener- 
ally more massive dark halo mass is needed for gas filling into 
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dark haloes and forming baryonic structure in the haloes. If 
T m = 7cmb the minimum halo mass for gas filling in is 
6.1 x 10 s Mq, while the minimum halo mass increases sig- 
nificantly for T m Tomb- In this sense, the formation of 
the first baryonic objects will be obviously suppressed in 
the small dark haloes located in host IGM gas. However, in 
the pure UCMH annihilation case without PBH radiation, 
we expect ~ 1O 6 M0 dark haloes attract baryons at z > 10 
in most cases expect for /ucmhw^qq > 10 . 

The annihilation energy deposited in a dark halo can be 
linearly divided into two parts: the energy from the back- 
ground where the cosmic UCMH annihilation occurs tb s d(z), 
and that within the local dark halo ei oc (z). The term eb g d(z) 
is obtained by equation (|39p . The local energy e\ oc (z) is a 
function of position inside the halo. We focus on the energy 
deposition at the center of the halo. The contribution by the 
extended dark matter in an isothermal halo is 

(av)c 2 (l 



Qo. 



2nm p m x f x 



fx) ( TTi \ 3 

CTtot^-C/xJPcorc 



T Rcov 



Ob 

-''•core 



where f x = IJdm/^m ~ 0.833. On the other hand, the lo- 
cal energy deposited by the UCMH annihilation depends on 
the UCMH distribution inside the halo. If we assume the 
UCMHs number density is uniformly distributed depending 



on the halo mass density, i.e., ar tYF I ' lI ^ ( \' const, (a relevant 



d "UCMH( r ) 

. rf/VfpMf.r) .. . 

distribution simulation see lSandick et al.ll201ll ). we have 
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[im p Mo M fx 



R~ 



3RI 



(49) 



In the following calculation we adopt equation (|49[1 for 
UCMH annihilation, and also include the extended dark 
matter annihilation within the halo. 

Figure [9] shows the gas evolution at the center of an 
10 6 Mq isothermal halo virializing at z v i r = 20 or z v - lr = 100. 
We also show the protohalo stage at z > z V ir- More energetic 
annihilation due to larger /ucmh or lower m x gives higher 
Xi on and /h 2 ■ The UCMH annihilation gives a significant im- 
pact on T m before virialization in the protohalo stage, that 
is because the cooling and heating mechanisms are different 
before and after virialization. The change of T m before viri- 
alization is mainly due the heating by background UCMH 
annihilation, which is totally dominated over the extended 
dark matter annihilation, thus brighter UCMH annihilation 
luminosity gives a higher gas temperature at z > z v j r . How- 
ever, after a dramatic temperature increase during the viri- 
alizing z ~ Zvir, H2 cooling becomes the main process to cool 
the denser gas at z < z v ir- The peak temperature during viri- 
alizing is around ~ 1000 - 2000 K. We find that higher H 2 
abundance, which is caused by brighter UCMH annihilation, 
gives a lower gas temperature after virialization for small z v i r 
(~ 20), but a higher temperat ure for large z v j r (~ 10 )- Thi s 
result is just between that in iBiermann fc Kusenkd (2006), 
who considered the effects of sterile neutrino decay can favor 
the structure formation, and lfUpamonti et alj l|200 7b). who 
showed that dark matter annihilation will slight delay the 
structure formation. The main reason of our difference from 
iRipamonti et all |2007bl ) for z v i r <C 100 is that, we take the 
baryon gas d ensity to be pr o portio nal to the halo density 
n b oc p as m iTegmark etafl (l997n . therefore more molec- 
ular gas due to stronger heating just means more efficient 
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Figure 9. Effects of UCMH and extended dark matter annihi- 
lation inside an 1O 6 M0 isothermal halo on the evolution of ion- 
ization (upper panels), temperature (middle panels) and H2 frac- 
tion /h 2 (lower panels) at the central region of the halo, where 
we choose the virial redshift z v ; r = 20 (thick black lines) and 100 
(thin blue lines), m x c 2 = 100 MeV (dash-dotted lines), 1 GeV 
(solid lines), 10 GeV (dashed lines) and 100 GeV (dotted lines). 



cooling. A more elaborate result can be made by adding 
more detailed gas dynamics and ener gy transfer including 
the UCMH radiation w ithin the halo (|Tegmark et al.lll997l ; 
IRi pamonti et al1l2007bl ). But such a new calculation should 
not change the fact that UCMH radiation, as well as the ex- 
tended dark halo annihilation, cannot change the gas tem- 
perature in the halo obviously after virialization. 

Note that the temperature T m in the halo only change 
by a factor of ~ 3 for a several orders of magnitude change to 
the UCMH annihilation luminosity inside the halo. There- 
fore, we cannot expect the first baryonic structure forma- 
tion can be obviously promoted or suppressed. After virial- 
ization, the effects of dark matter annihilation are always 
secondary compared to the H2 cooling mechanism. On the 
other hand, gas chemical properties such as X[ on and /h 2 
can be changed significantly that higher Xion and /h 2 are 
produced by brighter dark matter annihilation. 



4.2 Gas accretion onto PBHs and X-ray Emission 

If we consider gas accretion onto PBHs, and take the PBH 
fraction as j7/pbh > 10 -11 (10~ 12 ) for M pb h = 10" 6 M© 
(10 2 Afo), the first baryonic structure formation will be dif- 
ferent. Fig. 1 101 shows the minimal dark matter mass for the 
critical case T m = T v i r with different 77/pBH and the char- 
acteristic radiation Ex = 10 keV. Different from the dark 
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Figure 10. The minimal dark halo mass at T z 
radiation with ry/pBH = 10 — 7 



lines), 10 (dashed lines), 10" 



(dash-dotted lines), 10 
10 (dotted lines) and z 



for PBH 
" 8 (solid 

m = 100 



(dark lines), 50 (red lines) and 20 (blue lines). The characteristic 
PBH radiation is Ex = 10 keV. The shaded area is for the case 
that the minimal halo mass for gas collapsing becomes Mp>M < 
6.1 X 1O 3 M0 due to a lower ambient gas temperature compared 
to the CMB temperature. 



matter annihilation heating case (Fig. [8]), the minimal dark 
halo mass dramatically increase after z m . For z v j r = 20, the 
minimal dark haloes increase to > 10 6 A/ Q for?7/pBH > 10" 8 . 
Moreover, if we combine the annihilation before z m with the 
X-ray emission after z m , the minimal dark halo mass for 
T m = Tvi r can be even larger. 

Moreover, we expect the gas accretion onto PBHs in 
the dark halo environment above the critical mass in Fig. 
[10] will be slightly different from that in the ambient IGM, 
because the baryon gas is denser within a dark halo than 
the ambient IGM, which leads to a different accretion rate 
and baryon fraction inside UCMHs compared to the IGM- 
located-UCMHs. The accretion rate Mtjcmh and baryon 
fraction /(, inside a PBH host UCMH should be higher than 
those outside the halo. Therefore it is more difficult for X- 
rays from PBHs to pass through the host UCMH without ab- 



sorption. The critical redshift 



for X-rays escaping from 



the UCMH baryonic environment should be slightly delayed 
inside the halo than that in the background zj}f l0 < zj^ gd . 
It is possible that in a period of time that X-ray energy 
injection and deposition within a halo is mainly from the 
background even after virialization. 

The evolution of the baryonic structure inside a 1O 7 M0 
isothermal dark halo, as an example, is showed in Fig. 1111 
We consider two models: the number density of PBH host 
UCMHs being uniformly distributed per halo mass as men- 
tioned in Section Pl.ll also the uniformly distributed UCMHs 
per halo volume inside the halo as d " u - CMH oc const. Differ- 
ent UCMH distribution inside the halo would give differ- 
ent baryonic evolution. The local energy deposition for the 
UCMHs' uniform distribution per halo volume is written as 



eioc 



iaccO"tot£'XPcorG(l ~ fx) 
finipfx Vhalo 



He- 



Rtr 



(50) 
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Figure 11. Effects of gas accretion onto PBH inside a 10 7 M(; 



isothermal halo with z v i r 



: 20, 



bgd 



: 50, 



: 30 (left panel) 



and 2^ CMH = 15 (right panel). The characteristic Ex = 10 
keV with ij/pbh = 10 -7 (solid lines), tj/pbh = 10 -8 (dashed 
line, i)/pbh = 10~ 9 (dotted lines), ??/pbh = 10 — 10 (dash-dotted 
lines). We choose a more massive halo 1O 7 M0 because more mas- 
sive haloes are favors for the structure formation in this case. In 
the left panel we take d ?" CMH 
uniformly. 



uniformly and right panel gpSMH 



where L aC c is the total X-ray luminosity due to gas accretion 
onto PBHs inside the first baryonic object. The main results 
in Fig. [11] is very similar to those of dark matter annihila- 
tion in Fig. [5] For low virialization redshift z v i r = 20, the 
gas temperature is cooler for higher X-ray luminosity, but 
T m inside the halo only changes by a factor of 2 to 4 for a 
four orders of magnitude change in the X-ray luminosity. An 
more obvious change of chemical quantities (xi OI1 , /h 2 ) than 
the temperature change occurs for different X-ray luminos- 
ity within the halo. This means the effects of X-ray emission 
form PBHs on the gas evolution inside a halo is very small. 
PBHs which uniformly distribute per halo volume gives a 
slighter cooler gas within the halo, as well as lower Xi on and 
/h 2 than the uniformly distributed UCMHs per halo mass. 
That is because the latter model makes the average distance 
of UCMHs to be closer to the halo center, and gives more 
effects to change the gas properties at the center. 

In summary, UCMH radiation including both annihila- 
tion and PBH gas accretion enhances the baryon chemical 
quantities such as Xi on and /h 2 inside dark matter haloes 
which above the minimal halo mass for T m = T v i r , but the 
impact of UCMH radiation on the temperature of first bary- 
onic objects is small (by a factor of several), which shows the 
change of first baryonic structure formation due to UCMH 
radiation is less important than the H2 cooling and dark 
halo virialization time. However, the new chemical condi- 
tions provided by UCMH radiation can be more important 
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to affect the later gas collapse and first star formation after 
the first baryonic object formation, because more abundant 
H2 and electrons acting as the cooling agents can cool the 
gas more efficient during the gas collapse process, and pro- 
vide a lower fragmentation mass scale and first star mass 
(e.g JStacv fc Brommll2007r i. As we mainly focus on the first 
baryonic structure formation and evolution, the detailed cal- 
culation of first star formation due to the changed gas chem- 
ical components should be investigated more detailed in the 
future. 



5 DISCUSSION 

5.1 Status of UCMH Radiation in Reionization, 
Other Sources 

A variety of cosmological sources can reionize and heat the 
IGM at different redshifts before z ~ 6. So far we showed 
that UCMHs, even merely occupy a tiny fraction of total 
dark matter mass, provide a new gamma-ray background for 
gas heating and ionization. Also, the X-ray emission from the 
accreting PBHs could change the IGM gas evolution history 
dramatically after z m <C 1000, where the value of z m de- 
pends on the masses of PBH, host UCMHs and dark matter 
particles. Furthermore, we investigate that both dark matter 
annihilation and X-ray emission from UCMHs can dominate 
over the annihilation of extended dark matter halo. There- 
fore UCMHs are also an important energy source in dark 
matter haloes before the first star formation. In this section 
we briefly review all candidate energy sources during the 
Universe reionization era 10 < z < z oq . In particular, we 
emphasize the importance of UCMH radiation among all of 
these sources in different times. 

In this paper we focus on the heating and ionization 
processes after the last scattering epoch z ~ 1000, but some 
interesting effects can be produced by primordial energy 
sources at earlier time z > 1000. For example, cosmic gas 
heating and CMB spectral distortion at z ICC < z < z oq pro- 
duced by PBH gas accretion can be used to constrain the 
PBH abundance that / PBH < 10~ 8 f or M PBH > }0 3 Mq 
in the absence of UCMH annihilation (Ric otti et al.l 120081 ). 
However, as showed in Section \2. 21 UCMH annihilation can 
be more important than PBH gas accretion in the very ear- 
lier Universe even though UCMHs are just beginning to grow 
at that time. Hotter cosmic gas heated by dark matter an- 
nihilation suppresses the PBH gas accretion to becomes the 
dominated sources to distort CMB, and even changes the 
cosmic recombination process as showed in Fig. [2] Comp- 
ton y-parameter is likely to be used to constrain the UCMH 
abundance based on the annihilation scenario in future work. 

The influence of dark matter annihilation or decay at 
z < 1000 on the IGM during the reionization era has been 
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However, the annihilation power can be strongly increased 
by UCMHs in the early Universe, as what we have discussed 
in this paper. 

The next commonly suggested sources of ioniza- 
tion and heating is the first dark objects (dark haloes), 
which formed approximately at z < 100. Dark matter 
haloes enhanced the overall cosmic dark matter annihi- 
lation density due to the dark matter concentration in 
haloes illliev et al.ll2005l; lOda et al.ll2005l; ICiardi et al.ll2006l; 
Chuzhovll2008l; iMvers fc Nusser||2008l;lNataraian fc Schwarzl 



to co nstraint the cross section of the dark matter interac- 
tion JCirelli et al.ll2009l : ICalli et al.llioogl : ISlatver et al1l2009l : 
iKanzaki et alj 20101 ) . The basic assumption is that the dark 
matter distribution is smooth and homogenous at z > 100. 



20081; iBelikov fc Hooperl 120091 ; iNataraian fc Schwarzl 1200a 
Nataraian fc Schwarzll2010t ). The mass distribution of dark 
haloes varies from very low mass at (~ 10~ 6 Mp)) to high 
mass (~ 10 12 Mm) llGreen et alJliooi iDiemand et ai1l2005l ; 
iHooper et al.l l2007f ). depending on the d ifferent damping 
scale s due to different da rk matter models (|Abazaiian et al] 
l200ll; iBoehm et al.l|2005] ) as well as the mass-halo function 
l|Press fc Schechterl 19741 : ISheth fc Tormenll 19991 ). If UCMHs 
collapse with the homogenous dark matter together, the 
UCMH annihilation flux could still dominate over the to- 
tal annihilation flux within the dark matter haloes, at least 
in massive dark haloes with /ucmh-A/dm 3> rrih- However, 
remember that small dark matter haloes contribute to a 
significant part of the total annihilation rate after struc- 
ture formation. The profiles of the earth-mass dark matter 
haloes and the gamma -ray flux due to annih i lation have also 
been studied recently (|Diemand et al.ll2005l ; llshivama et alj 
2010). Similar to the UCMH emission, a large enhancement 
of annihilation signal is also expected due to the emission 
from the dark matter subhaloes as the remnant of structure 
formation at z < 60. Whether UCMHs or small haloes are 
more important for ionization and heat after z ~ 60 should 
be investigated in the future. 

The following ionization sources are the accreting 
PBHs, which locate in their host UCMHs, as mentioned in 
Section [23] PBHs with host UCMHs lead to a faster accre- 
tion than naked PBHs, but also absorb the X-ray emission 
due to baryon accumulation within the UCMH. The PBH 
accreting could only be more important than UCMH anni- 
hilation at z < z m < 1000 with sufficient abundance and 
radiation efficiency. Keep in mind that the X-ray emission 
here is from PBHs, or say, the PBH-UCMH systems, which 
paly an earlier role than the so-called accreting "first black 
holes (BHs)" , which are the remnants of first stars at 2 ~ 15. 

As mentioned in Section 2] dark matter annihilation 
or X-ray emission affects the baryonic structure formation 
and evolution. Also, they affect the process of first star for- 
mation. The standard fi r st star formation carried out at 
z ~ 20 l|Abel et al.ll2002l ; iBroom et al.|[2009h . but the first 
star formin g history can be affected by the primordial mag- 
netic fields l|Tashiro fc Sugivamall2006T ) , or by extended dark 
matter annihilation in the halo (so-call ed "dark star", see 
ISpolvar et all 120081 ; ISpolvar et alj|2009h . Previously it was 
said that the first stars gave the first light to end the cosmic 
"dark age", that cannot be true if exotic sources such as 
dark matter annihilation and accreting PBHs are included. 

Next ionization sources are more familiar to us. 
First stars emitted UV light and produced the "ion- 
ized bubbles", which could directly partially ionize the 
Universe at z < 20, or affected the coming forma- 
tion of ne xt generation stars and later galaxy f orma- 
tion (e.g., lHaiman fc Loebl 1 19971 ; IWvithe fc Loebl 120031 ; 
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IShull fc Venkatesaij|2008l ; IWhalen et~al1l2010l ). The death of 
first stars, produced the "first generation BHs", which emit- 
ted X-rays a nd ionized the Univ e rse at z ~ 15 or even closer 
dCenl 120031; iRicotti fc Ostrikerl 12004 iMadau et al l 12004 
iRipamontil 120071 ; iThomas fc Zaroubl 120081 ). Thereioniza- 
tion process was completed after galaxy formation, as galax- 
ies are generally considered the main candidates fo r the 
reionization of the Universe at z ~ 6 (Mciksin 2009|, and 
references therein). 

Future work that can be done includes studying the 
heating and ionization processes at 2 > 1000 due to an- 
nihilation, comparing the total annihilation rate from small 
dark haloes (1O _6 M0 < A/dm < /ucmh™' 1 as we mentioned 
above) with that from UCMHs, and distinguish the impacts 
of different ionization sources using the CMB polarization 
anisotropics and 21 cm spectra observational constraints. 
Actually, CMB polarization and hydrogen 21 cm line are 
powerful potential probes of the era of reionization to con- 
strain the early energy sources. The high multipoles of po- 
larization anisotropies may be able to distinguish UCMHs 
from small dark structures formed at z < 100, and further 
constrain the UCMH and PBH abundances. 



5.2 Different UCMH Profiles 



Remember that in Section [2.11 although several UCMH an- 
nihilation rate due to various profiles were given in Fig. [1] 
we choose the UCMH profile as p oc r -9 / 4 with a cut off 
at p m ax oc (t — ti) . Such a profile is based on the ana- 
lytical solution of t he radial infall onto a central overdensity 



I Berts chingcr 19851). A shallower density profile p(r 



which is given if the central accretor is a black ho l e, or a 
steeper profile p(r) oc r~ 3 simulated bv lMack et al.l l|2007l ). 
will change the total annihilation luminosity of a UCMH 
significantly. Fig. [12] gives an example of the different anni- 
hilation luminosity due to different density profiles in an en- 
tire UCMH. More concentrated dark matter distribution in 
a steeper profile leads to much higher total annihilation rate 
within the UCMH, because the center region of a UCMH 
contributes to most part of the total annihilation rate. How- 
ever, we conclude that the overall cosmic annihilation lumi- 
nosity density equation (I10p will not be changed too much 



for two reasons. First of all, the p 



profile usually 



appears in the outer region of a UCMH, but the change 
of the UCMH density profile at the outer region r 3> r cut 
will not dramatically change the total annihilation rate, the 
density distribution at the central region is crucial to deter- 
mine the total annihilation rate. Second, the contribution of 
PBH host UCMHs (with the profile p{r) oc r" 1 ' 5 near the 
center) to the overall cosmic annihilation should be much 
less important than the initial overdensity seeded UCMHs 



G°(r) 



both due to their shallower inner profile 



and the much lower abundance /pbh <S /ucmh- The anni- 
hilation luminosity should still be taken into account if the 
dark matter particle inner trajectory is high eccentric with a 
much closer pericenter than the cut off radius as in equation 
Q, but it is still lower than the luminosity of the overdensity 
seede d UCMHs as showed in Section 12.11 l|Lacki fc Beaconl 
l2010h . 




redshift 



Figure 12. An example of IGM evolution caused by UCMH an- 
nihilation with different UCMH profiles. In order to see the ef- 
fects more clearly, UCMH fraction is taken as /ucmh = 10~ 4 . 
The same group of lines from the top down show the results for 
m x c 2 = 1 GeV, 10 GeV and 100 GeV with p oc r" 9 / 4 (solid 
lines), p ~ r~ 3 (dashed lines) and p ~ r~ 1-6 (dotted lines). 



5.3 Extended Radiation Spectral Energy 
Distribution 

In Section [3] the products of two dark matter particles an- 
nihilation are assumed to be two gamma-ray photons both 
with energy m x c , and X-ray photons emitted from PBH 
host UCMHs are characterized with the single energy reflect- 
ing the center PBH mass E x 3 keV (M PB h/M©)~ 1/4 . The 
simplified treatment with a monochromatic (i.e., 5-function) 
spectrum of local UCMH emissivity is a good approach to 
demonstrate the crucial UCMH effects on the IGM evolution 
depending on the most crucial parameters such as UCMH 
and PBH abundances /ucmh and /pbh, as well as dark mat- 
ter particle mass m x . In this section we will experiment other 
extended photon spectral energy distributions (SEDs) for 
UCMH radiation rather than 5-function, and study the de- 
pendence of UCMH emissivity with SED. We will see that 
more elaborate considerations can quantitatively change the 
ionization results, but will not change the basic conclusions 
of UCMH radiation qualitatively. We first adopt a power 
law spectrum as F(E,z) oc E~ a , and then discuss another 
spectrum as F(E,z) oc E b exp(— dE), which are two most 
commonly used SEDs for annihilation and BH X-ray emis- 
sion. 

First of all we give an analytic calculation for UCMH 
radiation with locally monochromatic spectrum, based on an 
approximation that the optical depth described as equation 
(|38[1 can be neglected t < 1. This condition is applied to 
the spectrum E 1 (Ex) 3> 10 keV. Under the approximation 
r< 1 equation (137[1 in Section [3] can be simplified as 



dn 
dK 



EoEj Hoy 



l(z) 



(1+ 2) 3 / 2 ' 



(51) 



where the emitted photon number n and luminosity l(z) 
can be applied to both annihilation (n ann , Z ann (^)) or X-ray 
luminosity (nx, l&cc(z)), Eo for E x or Ex respectively, and 
l(z) oc (1 + zp with £ « 4 for annihilation and £ w 3 
for X-ray emission. We take the interaction cross section 
between photons and the IGM gas as atot (E) oc E~ k . The 
total energy deposition per second per volume equation (|40l) 
is integrated as 
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- C - k j l(z)n A (l+z) 3/2 a to t (-Bo), (52) 



«( 2 ) = 



where the subscript 5 is marked for monochromatic SED as 
in equation (|40|) . For high energy gamma-ray photons we 
can take k m 2 for Klein-Nishina cross section but for X-ray 
photons we take k ~ for Ex S> 1 keV. Another thing we 
mention is that ionization rate I(z) oc e(z) for £j on 1, so 
we only track e{z) based on different SEDs. 

The first extended SED is F(E, z) oc E' 1 for Ei < E < 
Ei, which gives a spectral number density for r ~ 



3/2 p -2 



(53) 



where A = E2/E1. Note that the number spectrum equation 
(|53[) is different from the equation (|51[) in the monochro- 
matic SED case. If we consider dark matter annihilation 
k ~ 2, the total energy deposition for UCMH annihilation 
is written as 



e(z) 



13 -2( 



fe In A 



1 + Z C 



1 + 2 



(54) 



Compared with equation (|52|) . the most obvious difference is 
the additional two factors (E /E 1 ) k and (1 + jz cq ) fe /(l + z) k 
significantly increase the annihilation energy density if Eo S> 
Ei. However This energy amplification for power law SED 
might be overestimated as we take the lower limit of the en- 
ergy (|40|l as E cq — Ei(l + z)/(l + z cq ), i.e, the Universe at 
redshift z can receive the emission from the matter-radiation 
equality era. A more general expression (Ei/Ei) can be 
used instead of (1 + z eq ) fc /(l + z) k with Ei being the thresh- 
old energy photons from UCMH radiation, then the am- 
plified factor compared to equation ()52|) can be written as 
(Eo/Ei) k 3> 1, no direct relation with E\ and E2. On the 
other hand, for the X-ray emission from a PBH host UCMH, 
as k « and the cross section is more or less the Thomson 
cross section, then the energy density is written as 



'-(*) 



l + ln(f)/lnA 



(55) 



In this case the energy density can be increased compared 
to the monochromatic spectrum if E\ > Ei, or decreased 
for Ei < Ei. But the logarithmic enhancement \n(Ei/Ei) in 
X-ray emission case is less significant than the annihilation 
case. 

For a more general power law spectrum F(E, z) oc E~ a 
(Ei < E < E2) with a > 1, the number density spectrum 
and energy deposition are calculated as 



dn 



dE~« HoVU^ \a+ ±± 



(!+*)-»/' W «(*) 



E-. 



EiE~ 



,(56) 



and 

e(z) 



l(z) 



HoV^m a + k — l \a + 



2 •» 
a+k—l 



71a(1 + z) 



3/2 



OC €s(z) X 



a+k-1 



(57) 



If we do not focus on the linear changed factors such as 

(a — l)/(a + -j C), mostly the power low spectrum will 

significantly increase the energy deposition e(z) as well as 
the ionization rate I(z) (and even a dramatic increase in 
some cases) by a factor of (E / E k ) k (Ei/ E i ) a+k ~ 1 . 

A black-body-like or say multicomponent spectral dis- 
tribution F(E, z) oc E b exp(-dE) with a peak Eo and c > 
can be approximately written as F(E, z) oc x h for x < 1 and 
F(E,z) oc exp(— dx) for x > 1, where x = E/Eo. In this 
case we find the number density spectrum as 



dn 
dEL 



cA 



l(z) (I + Z)- 3 / 2 




1) 



— — c 



(58) 



where A is an algebraic factor A — [(1 + b)^ 1 + d^ 1 ]^ 1 , and 
< f(d) < 1 can be obtained numerically, which is not im- 
portant for the following discussion. The terms in number 
density spectrum are proportional to E^ or E^ 1 ^ 2 . Accord- 
ing to equation (|58p. the integrated energy density e(z) is 
expected to be similar as equation (|52|) for b > k— 1. Other- 
wise the enhancement should be e(z) oc eg(z)(Eo/Ei) k ~ b ~ 1 . 
A steep x c with c > 1 for annihilation k ~ 2 and all the 
black-body-like SED for X-rays k ~ 0, are more or less sim- 
ilar to the 5-function SED at Eo. 

Now we study the case of Ex <SC 10 keV for X-ray emis- 
sion from very massive PBHs. Photon absorption in the IGM 
is important for Ex ~ 1 keV. The mean free path of X-ray 
photons describing by redshift change Az as 



T ~ Az- 



r 

or we have 

Az ~ 0.40(1 +z) 



-_(l+zr 5/2 n A (l+z) 3 a tot (lkeV) , 



-1/2 



Ex 
IkeV 



(60) 



where k ~ 3.3 for the photonionization cross section. Sim- 
ilarly as the former calculation, they energy density for a 
5-function spectrum is 



e s (z) 



0.40 



-l^ c (z)nA 



Eg 
IkeV 



(Ttot(-Bo) 



(61) 



Note that now we have a shallower density evolution e$ (z) oc 
(I+2) 3 compared with the transparently propagation for the 
case of equation (|54[) as eg(z) oc (1 + z) ( ' +3 ^ 2 . For a power 
law SED F(E,z) oc E~ a (Ei < E < E 2 ), we obtain the 
energy density as 



0.40c 




nA&tot(Eo) 



E 
IkeV 



(El 
\E 



(62) 



Therefore the energy density is enhanced for Ei > Ei, but 
the enhancement is less significant compared with that of 
equation (|57[) for the same spectrum index a. 

For a brief summary, the strength of UCMH radiation 
depends on its extended spectral energy distribution, which 
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will increase the properties of IGM ionization I(z) and heat- 
ing e(z). Compared to the basic results with the monochro- 
matic spectrum E x for annihilation or Ex for X-ray emis- 
sion, power law spectrum oc E~ a with a > can increase 
the energy density effectively, but black-body-like spectrum 
is more like the monochromatic SED case. For locally heat- 
ing Ex ~ 1 keV, the heating increases less significant than 
the transparently propagation case r< 1. Also, power law 
spectrum changes e(z) and I(z) more significant for annihi- 
lation than X-ray radiation. 



10* p — r-, 10- 2 




5.4 More Massive UCMHs Inside First Dark 
Haloes 

In Section [4] we assume that the number of UCMHs inside 
a dark halo is so huge that UCMHs in the halo are uniform 



distributed per halo mass 



= const, or per volume 



di\/ DM 

d"ucMH _ cons i This assumption can be invalid if the mass 
seed of a single UCMH raj is comparable to /ucmh-Mdm- 
In this case the position of each UCMH is important to 
determine the energy deposition within the halo. Usually 
X-ray emission can be neglected in a massive UCMH due 
to photon trapping (Section I2.2.3[) , we only focus on the 
UCMH annihilation. Equation (|49)l in Section I4.2l s invalid 
for massive UCMHs van ~ /ucmhMdm- In this case ei oc ,iso 
in equation (|49[) should be calculated as the summation of 
individual UCMHs 



£ioc,ucmh(z) — 



Q-tot(l - fx) 

fim p f x 



\ -> L V cMH,ip(ri) 
Z—i Airr 2 



(63) 



where Lucmh,; and r; are the annihilation luminosity and 
position of the i-th UCMH within the halo. 

An extreme case is that there is only one massive 
UCMH inside a dark halo, the energy deposited in the gas 
at the center of the halo can be obtained by equation (|63p 
in Section 14.11 Remember that for a uniformly distributed 
UCMHs with a same total mass of the single UCMH inside 
a halo we use equation (|49p . The ratio factor between (|63p 
and J49J is 



factor — 



B, c . 



ro 



3-Rtr 

RcOYC 



- 2 



Re 



tr 



3 (R c ore 

4£ V r o 



(64) 



where ro is the location of the single massive UCMH from 
the halo center. If we take R cola = £-Rtr with ^<1, the fac- 
tor can be approximately written as (3/4£)(7? core /ro) 4 . For a 
single UCMH located close to the isothermal case ro — -Rcore, 
the factor flM} is 3/(4£). For a single UCMH located at the 
turnaround radius ro — -Rtr equation (|64|l becomes ~ 
Taking a typical value £ ~ 0.1, the energy deposition ei oc 
contributed by a single UCMH varies from ~ 10 higher to 
~ 10 -3 lower compared to that contributed by the same to- 
tal mass but uniformly distributed small UCMHs. Fig. 1131 
gives the gas Xi on and /h 2 inside a 10 6 M© isothermal halo 
with a single UCMH located in difference locations. UCMH 
annihilation on the halo center will be less important for 
further located UCMH > 2R colc . The uniformly distributed 
UCMHs are more or less equivalent to a singe UCMH lo- 
cated at approximately ~ 1.5-R C orc- As an example in this 
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Figure 13. Ionization fraction x- lon (thick lines) and molecular 
hydrogen fraction H2 (thin lines) in a 10® Mq isothermal halo with 
a single massive UCMH located at -Rcore (dashed lines), 2i? C orc 
(dotted lines) and 10-Rcorc (dash-dotted lines) from the center of 
the halo. The solid lines shows the results given by many small 
UCMHs which number density is proportional to the halo density. 
Both single UCMH or small UCMHs is adopted a total UCMH 
fraction as /ucmh = 10 — 4 , z v i r = 100 and m x r? = 10 GeV. 



figure, we take the UCMH fraction as /ucmh = 10~ 4 and 
m x c 2 — 10 GeV. More massive dark matter particles or less 
UCMH fraction will give a faster decrease in energy depo- 
sition with the single UCMH moving outward the halo. We 
also check the importance of a single UCMH's position on 
the gas temperature inside a halo, but we see that the effect 
can be neglected for the temperature. 

On the other hand, the energy deposition by a single 
UCMH compared with the volume uniformly distributed 
small UCMHs with a same total mass can be amplified by 
a factor 



factor 



1 

6 

] 

6£ 3 



ro 
Be 



-Rtr 
Bcor 



1 - 



B c . 



Btr 



ro 



(65) 



which gives an enhancement factor about ~ l/(6£ 3 ) ~ l/£ 2 
for r = -Rcoro and a weaker factor f/6 for r = -R tr - These 
results show that the massive single UCMH is always more 
important for energy deposition compared with the volume 
uniformly distributed UCMHs. 

Therefore we conclude that different UCMH distribu- 
tion with a same total mass but different individual mass 
can change the energy deposition and structure evolution in- 
side a halo. More concentrated distribution towards the halo 
center or a closer located single halo near the center gives a 
more significant effect on the gas ionization and heating at 
the halo center. 



6 CONCLUSIONS 

Ultracompact Minihaloes (UCMHs) have been proposed as 
the primordial dark matter structures which formed by dark 
matter accreting onto initial overdensity or primordial black 
holes (PBHs) after matter-radiation equality z cq ~ 3100 
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l|Mack et alj|2007l ; iRicotti fc Gouldl [2009h . The key differ- 
ence between UCMHs and the first dark halo structures is 
that, UCMHs are seeded by primordial density perturba- 
tions produced in the very early Universe such as the phase 
transition epoches (10~ 3 < 5 < 0.3 for the initial overden- 
sity or 5 > 0.3 for PBHs), so they can grow shortly after 
2 ~ 2 cq . The radiation from UCMHs in the early Universe 
includes dark matter annihilation from all UCMHs, and X- 
ray emission from gas accretion onto PBHs. In this paper 
we investigate the influence of UCMH radiation on the early 
ionization and thermal history of the intergalactic medium, 
and the following evolution of the first massive baryonic ob- 
jects. Our conclusions are as follows. 

1. UCMH annihilation can totally dominate over the 
homogenous dark matter background annihilation, and pro- 
vide a new gamma-ray background even for a tiny UCMH 
fraction /ucmh = fi(z cq )/f2DM ~ 2.2 x 10 _is m~%o(l + z) 2 
with m Xi ioo = m x c 2 /100 GeV. We conclude that the influ- 
ence of dark matter annihilation on the IGM evolution can 
be significantly enhanced when we include UCMHs besides 
the homogenous dark matter background. In most cases 
the UCMH annihilation had been the dominated sources 
of ionization and heating gas since matter-radiation equal- 
ity epoch, until the X-ray emission from PBHs or large scale 
structure formation become important at z < 100. 

2. The impact of UCMH annihilation on the IGM can be 
approximately estimated by the quantity m~ 100 /ucmh- The 
threshold of UCMH abundance /ucmh to affect the IGM 
evolution by dark matter annihilation is m~\ 00 /ucMH > 
10~ 6 . After matter-radiation equality epoch, the IGM ion- 
ization fraction Xi on can be increased from Xi on ~ 10~ 4 in 
the absence of any energy injections to an upper bound 
Xion ~ 0.1, and the IGM temperature from the adiabati- 
cal cooling T m oc (1 + z) 2 to a maximum value T m ~ 5000 
K for the upper bound case m~\ 00 /ucMH ~ 10~ 2 , which 
is constrained using the CMB optical depth at late times 
z < 30. UCMH annihilation is able to significantly increase 
the Thomson optical depth r > 0.1 in the early Universe 
z 3> 30, which is unrelated with the measured CMB opti- 
cal depth at z < 30. The UCMH annihilation luminosity 
is based on the UCMH profiles, w here we take p oc r -2 ' 25 
from literature ( Bcrtschingcr 1985), steeper (shallower) pro- 
files decrease (increase) the allowed upper limit of /ucmh, 
but the variations of the overall IGM chemical and ther- 
mal quantities should not be changed too much, because 
the fraction of UCMHs with a profile p oc r~ as the PBH 
hosts are very small, and p oc r -3 occurs only in the halo 
outskirts r ^> r cut . 

3. Each PBH is located in its host UCMH (|Mack et all 

120071 ) . We emphasize that the impact of X-ray emission from 
PBH host UCMH systems is limited by the low abundance 
of PBHs (/pbh <fC /ucmh), the average inefficient radia- 
tion (rj <C 1), the photon trapping effect by the accreted 
baryons in host UCMHs, and outflows produced by rapid 
accretion feedback. Sufficient massive host UCMHs can ac- 
crete and thermalize the infalling baryons, which are ac- 
cumulated inside the UCMHs with a mass fraction of the 
UCMH f b > 10~ 3 , and trap X-rays from the accreting PBHs 
until a critical redshift z m ~ 32(Sm/MpBn) 1 ^ 2 i r n~ \qq , be- 
low which X-rays from a super-Eddington accretion flows 
onto PBHs could escape the surrounding baryon environ- 



ment in the host UCMHs. Although the PBH abundance 
is /pbh <C /ucmh due to the much higher perturbation 
threshold for the PBH formation, X-ray emission could dom- 
inate over UCMH annihilation and become more promis- 
ing cosmic energy source of the IGM ionization and heat- 
ing at z < z m if the PBH abundance is above a threshold 
»)/pbh//ucmh ~ 3.1 x 10~ 8 m^too(l + z ), which is only al- 
lowed beyond the standard Gaussian density perturbation 
scenario. 

4. As UCMHs are expected to exist in our Galaxy, we 
expect that UCMHs collapse with the homogenous dark 
matter background to form the first large scale dark mat- 
ter objects (dark haloes). If this is the case, the dark mat- 
ter annihilation from UCMHs inside the first dark halo still 
dominates over the extended dark matter annihilation back- 
ground inside the halo even after the halo virialization. 
UCMH radiation, including both dark matter annihilation 
and accretion emission, can dramatically suppress the for- 
mation of the low mass first baryonic structure, since UCMH 
radiation heats the IGM and provide a hot ambient gas envi- 
ronment up to T m ~ 10 4 K. The UCMH radiation enhances 
the baryon chemical quantities such as Xi on and /h 2 by or- 
ders of magnitude from Xion ~ 10 -6 and /h 2 ~ 10 -4 to the 
upper bound of Xi on ~ 10~ 4 and /h 2 ~ 5 x 10 -3 . However, 
the impact of UCMH radiation on the baryon temperature 
of the first baryonic objects is very small, which shows that 
the the influence of UCMH radiation on the temperature 
of first baryonic objects is small compared to the molecular 
hydrogen cooling and virialization time z v ir- However, the 
higher abundant Xi on and /h 2 provided by UCMH radiation 
decrease the gas temperature in the later gas collapse phase 
and can produce lower fragmentation mass scale and lower 
mass first stars. 

Also, we point out that, different spectral energy distri- 
butions of UCMH radiation also affect the processes of ion- 
izing radiation and heating gas. More concentrated UCMH 
distribution within a dark matter halo provides a more 
promising ionization phenomenon of the gas in first dark 
haloes. UCMHs should be distinguished from the small 
dark structure which formed during the structure formation 
epoch after z ~ 100. Future work need to be done to investi- 
gate the importance of small dark matter structure down to 
the earth mass compared with UCMH radiation in the early 
Universe at z ~ 60. Also, the CMB polarization anisotropics, 
21 cm spectrum and Compton y-parameter affected by the 
UCMH radiation also need to be further studied for a better 
constraint on the UCMH abundance. 
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APPENDIX A: PROFILES OF FIRST DARK 
STRUCTURE 



The dark 
as follows 
found 



m 



matter profile we use 
Most of the equa ti ons 



in this paper are 
list be lo w ca n be 
Tcgmark ct al. (|1997T ). iRipamonti I (|2007h and 
IRipamonti et al.l (|2007ar K A dark matter halo with a mass 
Mdm is assumed to distribute inside a truncation radius R tI , 
which is given by 



Rtr(z, Z v ir) = < 



f r 3 m dm i 1 / 3 

r 9 _ t( Z )-t(z ta ) i < < 

Vlr [ Z t(*Tlr)-*(*ta)J V ' r — — ta 



Rvir 

R\ 



(Al) 



z < z v 



Here z v ir is the redshift of halo virialization, zt& — 1.5(1 + 
z v ii) — 1 is the turnaround redshift. The dark matter mass 
inside the halo Alow is given as a parameter here. Dark 
matter within the halo is assumed to distribute uniform at 
z > Zta with a density evolution 



/9th (2, Zvir) = 

where A(z) 
given by 



pDM(z)e 



1.9A/(1-0.75A 2 ) 



(A2) 



(1 + 2 v ir)/(l + z). The virial radius Rvir is 
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3 M E 
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For an isothermal halo profile with a core radius J? C ore inside 
Rtr, the dark matter density p(z) is a constant p C orc for 
r < Rcotc, p(z) oc r~ 2 for -R CO rc < r < -Rtr- The density goes 
to the background dark matter density for r > Rtr. Here the 
core radius i? CO ro is obtained as 

R%i{z) Z > Zta 

Rcorc^Z , Zvir) — ^ Rvir (2 T^z )— £(zt yj Z v ir ^ ^ ^ Zta 



£i?vir 



•2 ^ Zvir 



where the coefficient £ is introduced as a parameter. And 
pcorc can be obtained integrating the halo mass within Rtr 
as Mdm- The total luminosity produced by the annihilation 
of dark matter particles distributed within the isothermal 
halo is 



J cxt,iso — 



2-Kpi 



(av)c 2 R 2 corc 



Rtr — — Rc 



(A5) 



On the other hand, the widely used NFW dark matter 
profile is p(r) oc r _1 (l + r-/7? corc )~ 1 for r < Rtr- Similarly we 
can write down the complete formula of density distribution 
as the function of redshift, and the luminosity produced by 
dark matter annihilation. 

The virial temperature of a halo T v i r is calculated as 



firtip GMom 

2k B Rvir '' 



(A6) 



sometimes people use the total mass Mhaio instead Mom/ x - 
However, as we consider if the gas from the background can 
collapse into the halo, we adopt the pure dark halo mass as 
to calculate the initial virial temperature. Therefore we can 



obtain the virial temperature as T v ; 



,2 x ltr 3 K (1 + 



2vir)(AfoM/M Q ) 2/3 , or Tvir = 380 K (Af DM /lO 4 M ) 2/3 (l + 
Zvir)/100. The baryon falling can occur when the IGM tem- 
perature Tm around the halo to be T m < T v i r , otherwise the 
gas pressure will impede the gas falling into the halo and 



collapsing to form smaller structure. Moreover, the condi- 
tion Tcmb = Tvir with Tcmb = 2.73 K (1 + Zvir) gives a 

" f . 



critical halo mass Mdm = 6.1 x 10 3 M<t 
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